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The California Digital Computer 


The California Digital Computer (CALDIC) is to be operated by the 
Engineering Department of the University of California at Berkeley for the 
benefit of its own research and that of the other departments and for the 
instruction of graduate students. This fact has strongly affected its logical 
design (in the hope of avoiding a large highly-trained staff) and also its 
physical design (requiring all parts to be accessible for demonstration and 
servicing). It is strictly a low-cost machine and makes no claims to superi- 
ority in speed or originality in design; most of the electronic techniques used 
are well established. It is not, however, designed for maximum simplicity 
of construction but rather for simplicity of description to and operation by 
the varied group of scientists who are expected to bring their computing 
to it. Its operation is to be straight-forward and as nearly like the familiar 
hand computational methods as possible, at least until experience is gained 
in its use. Its construction permits easy modification in the future, and it is 
expected that such modification will occur. 

The frontispiece shows most of the portions of the computer which are 
completed at this writing (September 1, 1950). At the left is the memory 
unit, consisting of a vertical rotating drum with four columns of combination 
reading-recording heads arranged around it and a set of electronic panels 
approximately six feet high by ten feet wide which contain the circuitry 
required to write into and read from the drum. To the right of these is the 
first of the arithmetic registers; the second has been built but not yet in- 
stalled on the blank panels below. At the extreme right is the order counter 
which is a shifting and counting register designed to store and supply the 
address of the next order to be carried out. The panels shown, a few of which 
are temporary, were built first to allow their operation as a test system, and, 
as other parts of the computer are completed, they will be added to the 
system. 

Figure 1 shows a simplified block diagram of the computer. It is a serial 
machine so far as decimal numbers are concerned; but the binary parts of 
each digit are transmitted in parallel, and most of the circuits have four 
channels. The three arithmetic registers are of the shifting type, with circu- 
lation paths provided as shown. The “A”’ or accumulator register incor- 
porates an adder in its first column, and circulation of its contents simul- 
taneously with receipt of a new number causes the sum to appear and stand 
in the register. For multiplication the ‘“‘D’’ register holds the multiplicand, 
which is repeatedly circulated but simultaneously shifted and added into the 
A register as required by the digits of the multiplier which are drawn from 
the “R” register. As the multiplier digits are used the overflow from the A 
register appears in their places, and at the end of the process the double- 
length product occupies both the A and the R registers. 

All timing pulses for the computer are drawn from a timing track perma- 
nently recorded on the drum, and therefore no effort to control the drum 
speed is needed. A small gap in the timing track is used to provide an “origin 
pulse,”” and recorded words are located by counting the timing pulses, be- 
ginning the count at the origin. Information recorded on the drum is there- 
fore not lost during power failures or other shut-downs; subroutines could 
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be permanently stored if desired. The timing pulses are divided by suitable 
circuits into “digit pulses’ in groups of eleven which provide for shifting 
numbers through the registers and “‘space pulses’’ which occur between the 
groups of digit pulses and which initiate all switching. 


COMPUTER BLOCK DIAGRAM 
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FIGURE 1 


The operation of the computer will be cyclic, essentially as follows: 


1. The address of the next order to be performed is shifted from the 
order counter to the address register and used to extract the order from the 
memory into the D register. 

2. The order is shifted into the operation and address registers, and 
simultaneously the contents of the address register are returned to the order 
counter and increased by one to provide for the next cycle. 

3. The address register now contains the address of the operand, which 
is withdrawn from the memory into the D register and operated on as re- 
quired by the contents of the operation register. At the end of the cycle the 
result stands in the A register. 


The normal cycle requires two trips to the memory, each of which might 
take almost the time of a full drum revolution (17 milliseconds) in addition 
to the time required for the operation itself which is usually small in com- 
parison. In most cases, however, the orders will be in consecutive memory 
boxes and the whole cycle will require only slightly more than the time of 
one revolution. In addition, the consecutively-numbered boxes are, spaced 
around the drum at intervals which correspond to the time required for an 
addition and order-shift, and recognition of this fact in programming may 
make it possible to carry out as many as ten operations per revolution of the 
drum. Abnormal cycles are used for several orders for which no new operand 
is needed; shift-left, shift-right, and sub-program orders are of this type. 
The cycle time for these, however, is about the same as for normal cycles. 
Multiplication and division, which are performed by repeated addition, may 
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require an additional revolution. Because of these variations the average 
operation speed is difficult to estimate, but it should be considerably better 
than 50 single-address operations per second on most problems. 

The orders to be provided are: add, subtract, multiply, multiply and 
round off, divide, divide and round off, extract the square root, read the 
input tape, transfer to memory, print out, clear address and add, change 
contents of order counter, change contents of order counter but only if the 
previous operation overflowed, shift left, and shift right. 

The only unusual orders are the “clear and add” and the conditional 
sub-program orders. The former provides for the modification of orders by 
building the desired address in the A register and attaching the operation 
part of the order to it. The latter provides for discrimination and also for 
overflow of the accumulator; a possible overflow is provided for by inserting 


INPUT— ELECTRIC TYPEWRITER 


a NUMBERS TYPED LEFT TO RIGHT 
(MOST SIGNIFICANT DIGIT FIRST) 


<+— PUNCH 
<—— PUNCHED PAPER TAPE 
*— PHOTOELECTRIC READER 


ORDER OF DIGITS REVERSED 


*—D REGISTER 


lo LEAST SIGNIFICANT DIGIT TRANSMITTED FIRST 


*—A REGISTER 








<—MAGNETIC DRUM 
FiGurRE 2 


the instruction which tells the computer what to do if the overflow occurs. 
If no overflow occurs, the instruction is disregarded. 

Orders consist of two decimal digits which specify the operation and four 
additional digits which usually specify the location in the memory of the 
operand, but sometimes they give other information such as the location of 
a sub-routine or the number of places of shifting in a shift-left order. No 
use is made of the sign column or the remaining four digit positions of the 
normal 10-decimal-digit word. 

Input will be from perforated tape, using the system diagrammed in 
Figure 2. In general all orders and data will be stored in the memory at the 
beginning of a problem, and the input will not be consulted further during 
the computation. The ‘‘read the input”’ order is designed mainly to instruct 
the machine to proceed to the next problem when one is completed; the tape 
reader will be loaded manually, and the tape will be run clear through on 
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receipt of the “‘read’’ order. The information read from the tape may, how- 
ever, be placed in any desired portions of the memory, as specified by ad- 
dresses punched into the tape ahead of the blocks of information. The entire 
memory can be filled in approximately three minutes. Output also will be 
through punched tape followed by printing with a standard tape-reading 
typewriter. Output speed is limited to punch speed of about ten digits per 
second. 

The outstanding feature of the computer is its magnetic memory, shown 
in more detail in Figure 3. The drum is about eight inches in diameter and 





Ficure 3. Magnetic Memory 


26 inches long and will store 10,000 ten-digit numbers with their signs. Since 
the no-excess binary code is used, over 480,000 memory cells are required 
at a density of about 900 per square inch of surface. There are 200 specially- 
designed magnetic reading-recording heads, and each is provided with a 
writing tube, a writing gate, and a reading gate. In addition 50 tubes are 
used in the band switch, and with the timing circuits and the 4-channel 
reading and writing amplifier circuits the memory unit requires about 750 
tube envelopes, many of them double triodes. This is more than half of the 
approximately 1300 tubes used in the whole computer along with 1000 
crystal diodes. 

Practically all of the design and construction of the CALDIC is the work 
of graduate and undergraduate students in Electrical Engineering, some 35 
having been associated with the project already. The electronic work is now 
being supplemented by research and instruction on programming and appli- 
cations. The need for both types of training is demonstrated by the demand 
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which has developed for graduates of the program. Financial support from 
the Mathematics Branch of the Office of Naval Research has made the 
work possible, and techniques developed by other ONR-supported projects 
have been used freely. 

Paut L. Morton 
Computer Laboratory 
University of California 
Berkeley 4, California 


Step-by-Step Integration of X = f(x, y; Z, t) 
without a “Corrector” 


Introduction.—For step-by-step numerical integration of ordinary differ- 
ential equations there are too many formulae, too few evaluations or com- 
parisons. Perhaps the complexity of the subject will not permit the general- 
izations the mathematician would prefer, but only restricted numerical 
comparisons of specific procedures. This paper, at any rate, will restrict 
itself to showing reasons for preferring one of two procedures for the integra- 
tion of second-order differential equations in which the second derivative, 
i= a , is a function of z alone, or of z and ¢, or of z, y, 2, ¢ (with similar 
equations for 7 and 2). These equations are important, of course, in certain 
dynamical problems of astronomy, ballistics, aerodynamics, etc., including 
the rocket problem. 

The comparison of the two procedures will extend to at least three 
equivalent forms, in which the integration formula at each step is based upon 


(a) antecedent values of 2, 
(b) backwards differences of Z, 
(c) central differences (estimated) of Z. 


Although there are reasons for preferring forms (a) or (b) in certain circum- 
stances, the comparison will be made first between the procedures in form 
(c) because there is some advantage to distinguishing between errors of 
estimation and errors resulting from the neglect of higher order terms in 
the integration formulae. 

The preferred procedure, designated the “‘second-sum procedure’ or 
“>? procedure,” involves tables of and its first and second sums, 2% and 
>*#, as well as its differences (explicitly or implicitly), 6%, #%, #®Z%, ---. It 
has been used for more than a century by astronomers,” but has apparently 
been overlooked by a number of mathematicians, physicists, and others, in 
recent years. 

The compared procedure, designated the ‘‘second-difference procedure” 
or “‘8? procedure,” involves the summation (explicit or implicit) of the second 
difference of x, &x, which is obtained by formula from ¢ and its differences 
(or the equivalent). This procedure was used by COwELL & CROMMELIN, 
perhaps for the first time in an extended calculation, in their celebrated 
prediction of the return of Halley’s Comet in 1910. In publishing the results 
of this work* they recommended without explanation, however, that future 
integrations of this type should be done with the 2* procedure and formulae. 
In spite of this recommendation a great deal of work has been done by the 
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& procedure, by astronomers as well as others, that could have been shown 
to be better adapted to the =? procedure. 

In short this paper will show that (1) under comparable circumstances, 
the & procedure requires two approximations (or the use of “predictor” and 
“corrector” formulae*) whereas the 2? procedure requires only one. It will 
show further that (2) the estimated differences normally employed by the 
astronomer may be replaced by backwards differences or antecedent values 
of the second difference without increasing the error or requiring a second 
approximation. The second of these facts, as well as the first, is significant 
to modern calculation with automatic electronic digital computers. 

The 8? Procedure.—For the purposes of the comparison it will be neces- 
sary first to outline the “8? procedure” with estimated central differences, 
since these may be unfamiliar to many persons otherwise well acquainted 
with the problem. Their introduction makes it possible to use the same for- 
mula for “predictor” and “corrector,” the difference between the two ap- 
proximations being due to the improvement in the estimates of the differ- 
ences involved. This formula’ we shall call the ‘‘(5.2) formula” and write 
as follows: 


2 =— J? a — _ — £42 om S68 .. oo 
(6.2) Fr, = h {n+ 20%, 740 ote + 50480 ° 2" 


where h is the interval of the argument (0.1 in the table below), the notation 
for the differences is best appreciated by reference to the numerical table 
below, and the subscript indicates that the various quantities are to be 
found on the same line of the table (the underlined quantities on the 0.8 
line in the example below). 

The following table shows an integration table for the simple differential 
equation ¢ = —z. How the table was started is not significant; at the given 
step all numbers not in parentheses may be taken as correct. The values of 
the differences of @ in parentheses are estimates based upon the assumption 
that the sixth difference is zero. The estimates of 59.3, 52.35, and Zo.9, on 
the other hand, are a result of the first approximation or “prediction” below 
the table. 





z bz 8x #=-2 bf aE ee se OF 
+.000 000 000 000000 —.0000000 0000 000 

+99 833 417 —99 8334 +9975 —100 
+.099 833 417 —997 502 —.099 8334 +9975 —100 

+98 835915 —98 8359 +9875 -91 
+.198 669 332 —1985038 —.198 6693 +19850 —197 

+96 850877 —96 8509 +9678 —100 
+.295 520 209 —2952740 —.295 5202 +29528 —297 

+93 898 137 —93 8981 +9381 -91 
+.389 418 346 —3890939 —.389 4183 +38909 —388 

+90 007 198 —90 0072 +8993 —89 
+.479 425 544 —4790 261 —.479 4255 +47902 —477 

+85 216937 —85 2170 +8516 -39 
+.564 642 481 —5641721 —.5646425 +56418 —566 

+79 575 216 —79 $752 +7950 (—89) 
+.644 217 697 —6436810 —.6442177 +64368 (—655) 

+73 138 406 —73 1384 (+7295) (-89) 
+.717 356 103 (—7167 586) —.717 3561 (+71663) (—744) 

(+65 970 820) 


(+.783 326 923) 
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“Prediction” “Correction”’ 
+ Zo.8 = —0.717 3561 —0.717 3561 
1 
+ 12 °es = +5971 .9 +5973 0 
ail 50.8 = +3 .1 +3 0 
240 0.8 = . ° 
h-*5*Zo.8 = —0.716 7586 .0 —0.716 7585 .0 
820.8 = —.007167586 — .007167585 


After the predicted 5*z9.3 is summed in the table to produce an estimate 
of 20.9, it is possible to calculate accurately #9... = —0.783 3269. [From an 
inspection of the “‘correction” we may verify that no significant revision of 
this number will be necessary. ] The table may now be differenced to show 
that the accurate value of &#o., = +71676, and this number replaces the 
estimate (+ 71663) in the “correction,” with a resulting change of one digit 
in 829.8, 50.85, and 2,9. The table is now ready for the next step. 

The =? Procedure.—The ‘‘(2,”) formula”’ is as follows: 





(22) wa = Wt (2, + , — aia Pe t sz, — | 


12 60480 


where the notation is the same as in the (4) formula. The numerical co- 
efficients are the same in the two formulae because the one is simply the 
second sum or the second difference of the other. The table associated in this 
procedure with the numerical example of the preceding section is as follows: 


t zz zz = -2z bz Fz eE 5 &z 

0.0 +0.000 0000 — .000 0000 0000 000 

+9.991 6653 —99 8334 +9975 —100 
0.1 +9.991 6653 — .099 8334 +9975 —100 

+9.891 8319 —98 8359 +9875 —97 
0.2 +19.883 4972 —.198 6693 +19850 —197 

+9.693 1626 —96 8509 +9678 —100 
0.3 +29.576 6598 —.295 5202 +29528 —297 

+9.397 6424 —93 8981 +9381 —91 
0.4 +38.974 3022 — .389 4183 +38909 —388 

+9.008 2241 —90 0072 +8993 —89 
0.5 +47.982 5263 — .479 4255 +47902 —477 

+8.528 7986 —85 2170 +8516 —89 
0.6 +56.511 3249 — .564 6425 +56418 —566 

+7.964 1561 —79 5752 +7950 (0) 
0.7 +64.475 4810 —.644 2177 +64368 (—566) 

+7.319 9384 —73 1384 (+7384) (0) 
0.8 +71.795 4194 —.717 3561 (+71752) (—566) 

+6.602 5823 (—65 9632) +6818 
0.9 +78.398 0017 (—.783 3193) (+78570) 


NoTE 1: The estimates of the differences (in parentheses) are arbitrarily 
made poorer in this example than in that of the preceding section, to empha- 
size the fact that an even larger error of estimation will not require a second 
approximation in the 2? procedure. 

Norte 2: (For those unfamiliar with sum columns.) Given a starting 
value, the numbers in the 2% column are summed from the ¢ column just 
as the # column could be summed from the 6¢ column. Similarly the 2*z 
column is summed from the 2Z column. Starting values may be obtained by 
an inversion of the process followed below. 
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Note 3: The subscript x refers to the 0.9 line in this example rather than 
the 0.8 line, as indicated by the underlining, in order to make the figures of 
the two examples more directly comparable. 

The calculation of 29.9 and of an accurate value of Zo.) to replace the 
estimate of this quantity proceeds as follows, in accordance with equa- 
tion (2,2): 


+ Bog = + 78.398 00 2 

+ i‘ ian Oo Se 

= oe -3 3 
7a Pts 





h? 20.9 = + 78.332 69 2 
to.9 = —Z%o.2 = — 0.783 3269 


Although this value of £o., differs by 76 units of the last place from the 
estimate in the table, it will be found, in checking the figures of a possible 
“correction” or recalculation, that the preliminary value of #;%o.9, and hence 
the final value of #o.9, will not be altered from what is shown in the above 
calculation. Nor would they if an even greater error were made in the 
estimate. 

Comparison.—1. It will be observed, in the numerical work of the two 
examples, that the #; term involves in the first a 5 digit number 
(#%o.3 = +71663), in the second a 7 digit number (Z%.9 = — 7833193), but 
that only 5 digits are carried in either calculation (+59719 and —65277). 
That is in the =? procedure the last two digits of the estimate of # are of no 
significance, and a much larger estimation error can be tolerated than in the 
8z that plays the same role in the & procedure. Thus the first approximation 
is sufficient in the 2? procedure, but a second is necessary in the & procedure 
to alter or check the last digit in the calculation. It is assumed, it will be 
noted, that successive columns of sums, function, and differences taper off 
in size by about one digit or more per column. When the tapering is more, 
the situation is even more favorable to the =* procedure: for example, if 
there are four more digits in the 2*¢ column than in the # column, the last 
four digits of the estimate of # will not be significant, and a very crude 
estimate is tolerable. A tapering of one digit per column is probably the least 
that would be tolerated in any procedure for numerical integration when 
automatic or rote methods are employed; an experienced computer, using 
the utmost of judgment, may be able to integrate by these procedures when 
the tapering is only one digit per two columns. 

2. Since a larger estimation error can be tolerated in the =? procedure 
than in the & procedure, it is evident that fewer difference columns will need 
to be taken into account, an advantage when the memory of an automatic 
computer is limited. 

3. It is assumed that the accumulation of rounding error and the effect 
of neglected terms will be about the same in the two procedures. The discus- 
sion is concerned instead with the effect of errors of estimation, explicitly 
as in the foregoing example or implicitly in equivalent processes involving 
backwards differences or antecedent values of #. (The estimation of differ- 
ences should not be confused with a first approximation or “‘prediction,” 
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but should be recognized as merely an alternative way of making the inte- 
gration depend upon backwards differences.) 

4. It is assumed that the purpose of numerical integration is the con- 
struction of as accurate as possible a table from which the integral may be 
obtained, and not necessarily a table of the integral itself. Often, as in comet 
orbits or ballistics tables, only a few values of the integral at the end of the 
table are needed. Since, however, it is customary to carry two or three more 
places in a table than are required in the end, it will be noted that the values 
of the integral obtained in the first approximation by either procedure are 
probably sufficiently accurate for all intended uses. From this point of view, 
evidently the second approximation in the & procedure is necessary only to 
avoid the accumulation of error in the summation of 6z. The equivalent 
summation of ~ in the 2* procedure is free of estimation error. For inter- 
polating non-tabular values of the integral, Bower has supplied a table.* 

5. To shorten the paper, I have not included in the foregoing discussion 
a third procedure which may bé used in integration of second-order differ- 
ential equations. This procedure is based upon the formula for the second 
derivative, #, in terms of the differences of x, which we may write 


1 1 
Gr = ht + 75 Ht — 95 Ot + . 


My investigation of this procedure revealed only a greater disadvantage: 
namely, three approximations as compared with the two of the & procedure 
and the one of the 2? procedure. Various attempts have been made to im- 
prove this formula for restricted classes of second-order differential equations 
by introducing special devices to eliminate the 7; term—beginning with 
NuMEROV,‘ and most recently by Fox & Goopwin.° At best these devices 
yield results comparable to the =* procedure only when the 1/240 term is 
negligible, and at least for some of them JACKSON’ is right in believing that 
they would become the 2? procedure if they were carried to their logical 
conclusion. 

Another proposal sometimes made in connection with the second- 
derivative formula above is that it be used to correct a first approximation 
in which the estimation error is ignored. This proposal would have merit if 
it were not for the fact that the estimation error is negligible in the first 
approximation by the 2 procedure. 

The Backwards-Difference Formula.—When the central differences used 
in the (2.2) formula are estimated by a summation based upon the repetition 
of the last value of a given difference, the calculation will yield exactly the 
same result as if it had been based upon a backwards-difference formula 
ending with the same difference. This formula is as follows: 


(2,?) Z, = h? {22 + 5 [a + bZn-3/2 +> ° 5! Fn—2 


18 1726 : 
20 Pte—w2 + 9016 Os + ig saa gia + * | - 


The estimation of differences may thus be entirely avoided in rote calcu- 
lation, by automatic machinery or otherwise. 


+= 








66 STEP-BY-STEP INTEGRATION OF ¢ = f(x, y, 2, #) 


The Antecedent-Function Formula.—When the computer feels that he 
can dispense with the check afforded by the run of the differences, he may 
replace the (2,?) formula by the equivalent expression in terms of ante- 
cedent values of the function to be eiioes 


(27) t=}? 22. +3 + : 5(1 | +> > + 5 4. ++) fa 


“Shae + 3- 1 


-) 
+(3 + 3-35 + +) aes 
) 


1 / 18 
~1 (384... 


It is evident that this general expression will yield a number of formulae, 
depending upon the order of the last difference of the (2,”) formula that is 
taken into account: 


(Za) zt, = h? {22 + r Ln-1 
(222)1 %,=h* {22 +7 + DR (2%,-1 — 29} 
(Za")2 Ln = h? {2, + m6 (59%,-1 - 58%n—2 + 192,..)| 


(Z4”)s Ly = h? {220 + 33 + —— 40 (17 ¥n-1 — 112%,-2 + 73%,_3 — 182,..)| 


and so on, where (2,”)» includes the effect of 52. 

I am greatly indebted to Professor W. E. MILNE and Dr. GERTRUDE 
BLANCcH, of the Institute for Numerical Analysis, for advice and comment 
on this paper. A part of the work was done while I was on the Institute staff 
and so had support from the Office of Naval Research. 

SAMUEL HERRICK 
Department of Astronomy 
University of California 
Los Angeles, Calif. 


1T. R. vON OPPOLZER, Lehrbuch zur Bahnbestimmung der Kometen und Planeten, v. 2, 
Leipzig, 1880. 
. C. Watson, Theoretical Astronomy. Philadelphia, 1892. 
=F CowELL, & A. C. D. CROMMELIN, Investigation of the Motion of Halley's Comet 
from 1759 to 1910. Appendix to Gr ich Observations for 1909. Edinburgh, 1910. p. 84. 
4B. V. NumERov, Méthode nouvelle de la détermination des orbites et le calcul des v7 hémé- 
rides en tenant c —_— des perturbations. Observatoire Astrophysique Central de Russie, 
Publications, v. 2, Moscow, 1923. “‘A method of extrapolation of perturbations,” Roy. Astr. 
Soc., Monthly Notices, v. 84, 1924, p. 592-601. 


5 J. Jackson, Note on the numerical integration of Fa pi = f(x, t). Roy. Astr. Soc., Monthly 
Notices, v. 84, 1924, p. 602-606. 
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6° E. C. Bower, “Coefficients for interpolating a function directly from a table of double 
integration.” Lick Observatory, Bull., no. 445, v. 16, 1932, p. 42. 

7 Interpolation and Allied Tables, —= from the British Nautical Almanac for 1937, 
4th edition with additions. (London, H. M. Stationery Office, 1946.) 

8W. E. Mine, Numerical Calculus. Princeton, 1949. 

*L. Fox & E. T. Goopwin, “Some new methods for the numerical integration of ordi- 
nary differential equations,” MS. awaiting publication. 


Formulas for Calculating the Error Function 
of a Complex Variable 


1. Various methods have been suggested for the computation, to high 
accuracy, of the error function #(Z) = /,%e-“du for complex arguments 
Z = X + iY. The Maclaurin series is convenient for extreme accuracy only 
when |Z| is small, and the asymptotic expansions' are useful only for fairly 
large values of |Z|. Other less elementary methods (e.g., the AtrEy con- 
verging factor,? or the continued fraction expansion’) have limited success 
in certain regions of the Z-plane. The most convenient methods for calcu- 
lating (Z) to very many figures are described in the works of MILLER & 
Gorpon? and Rosser.‘ It is the purpose of this note, which is self-contained, 
to present in a concise and practical form, two schemes for the computation 
of &(Z) which follow the main ideas of Miller, Gordon, and Rosser. 

2. The basis of the present methods is the following formula: 


(1) >< exp (—(u + na)’?) = ria S exp (—n*x*a-*) cos (2n2u/a) 
which is an immediate corollary of Poisson’s formula.* Formula (1) is also 
essentially Jacobi’s imaginary transformation,® which has long been familiar 
to a number of mathematicians and physicists.’* Also, (1) was used by 
Dawson’? in his computation of fo” e“du. Two different methods for calcu- 
lating @(Z) are given here in 4 and 5. In all formulas, summations are from 
1 to © except where otherwise noted. The symbol = will denote approximate 
equality. 

3. Dawson’s formula for /o¥ e“du is needed in the second method. It 
follows from the approximate equality obtained from (1), which is 


(2) e@(1 + E) = ax-*(1 + 2 ¥ exp (— a’n’) cosh 2nau), agi 


where the relative error E = 2 > e~”**/ cos (2mxu/a), is of the order of 
magnitude of 2e-**/**, since 2e~***/* is very small by comparison. Approxi- 
mate values of E are given in the following table: 


a 1 0.9 0.8 0.7 0.6 0.5 
E 10~ 10-° 3-10-* 3-10-* 3-10 4-10-" 
Putting a = } in (2) Dawson obtains 
(3) e* = w-4(4 + Y e-”* cosh nu), 
the relative error in (3) being less than 2-10-'’. Then integration of (3) 


results in 


(4) f 3 e“du = x-4(4Y + & mn exp (— n*/4) sinh nV). 
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4. One method of obtaining a formula for #(Z) is to choose the path 
of integration from 0 to Z to consist of the line segments [0, X] and 
[X, X + 7Y]. We find 


x 4 
(5) #(Z) = f e~“du + anf e“ sin 2Xudu 
0 0 
4 
+ ie-** f e“ cos 2Xudu. 
0 


The first term on the right is numerically well known.'® To deal with the 
other terms, we multiply across in (3) by sin 2Xu and cos 2Xu respectively, 
and integrate between 0 and Y. Using the explicit expressions for integrals of 


the form f cosh nu po : du we readily obtain (Z) = A + iB, where 


x 
(6) A= f e~“du + x-te-**(1 — cos 2X Y)/4X 
0 


+ ate SY (n? + 4X*)—"(2X — 2X cosh nY cos 2X VY 
+ nsinh n¥Y sin 2X Y)-exp (— n?/4) 
and 


(7) B= w-te-**((4X)— sin 2X Y + DY (mn? + 4X?) exp (— n*/4) 
- {2X cosh nY sin 2X Y + nsinh nY cos 2X Y}). 


5. Another formula for calculating 6(Z) is found by choosing the path 
of integration to consist of the line segments [0,iY] and [7Y, X + iY], 
whenfwe find 


Y x 
(8) #(Z) = if e“du + e¥? f e~“ cos 2 Yudu 
0 0 
xX 
— ie¥* £ e~ sin 2 Vudu. 
0 


The first term on the right can be found from Dawson’s approximation in (4). 
The other two terms are found by letting a = 12 in (1), then multiplying 
across by cos 2 Yu and sin 2 Yu, and integrating between 0 and X. In view 
of the fact that /° e~“du is negligible, being of order 10—”, it is easily shown 
that if we neglect the integral of everything arising from m ¥ 0 in the left 
member of (1), when X < 6 the error is < 2 {4% e-“du. When X > 6, we 
can replace X by 6 in (8), with a relative error < /;* e~—“du, and in view of 
the result in the preceding sentence, the total error is surely < 3 (4° e—“du. 
For maximum accuracy it is not necessary to take more than 24 terms 
arising from under the summation sign in the right member of (1), 
because then the truncating error, which is of the order of e~**’, is less 
than the error in the approximation formula. Making use of the explicit 


expressions for integrals of the form f cos pu a du, whenever 
| Y| ¥ nw/12, we obtain #(Z) = C + 1D, where 
Pad 


24 
e¥*((2Y)— sin 2X Y + 72 } (m*x? — 144 ¥*)— exp (— n*x*/144) 


n=1 


{gan sin $xnX cos 2X Y — 2Y cos §rnX sin 2X Y}) 


9) cae 
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and 
(10) D =r (3Y + ¥ mn exp (— n’*/4) sinh nY) 
1 “4 
-- mie?” (2 Y)-"(1 — cos 2X Y) + 72 ¥ (m*x* — 144Y?)- 


n=1 


X exp (— n*x?/144) | = sin inwX sin 2X Y 
+ 2Y cos ixnX cos 2X Y — 2 r}) 


In the exceptional case when | Y| = nx/12 for some value of n, omit that 


| 


term from inside the summation signs in C and D and instead add 
(11) mwi(48VY)-'(4X VY + sin 4X Y — i{1 — cos4XY}). 


Both methods which are described in 4 and 5 are accurate to within a relative 
error of about 10-'* in |@(Z)|, éxcept for values of | Y| very close to m/12 
in the latter method. The amount of computational work seems to be about 
the same for these two methods, and in general it appears difficult to recom- 
mend one of the pairs of formulas (6), (7) or (9), (10) in preference to the 
other. Each method affords a good computational check upon the use of 
the other. 

6. Often the value of e?*@(Z) is desired because it is a smoother function 
which can be tabulated to a constant number of decimals. Thus from (6) 
and (7) we have at once e?°@(Z) = E + iF, where 


x 
(12) E = exp (X? — Y?) cos oxy f e~"'du 
0 


+ a-e-¥*((4X)—"(cos 2X Y — 1) 


+ 2X & (n? + 4X") exp (— n?/4) {cos 2X Y — cosh nY}) 
and 


x 
(13) F= exp (X? — Y’) sin oxy f e~“du + w—te-**{ (4X)—' sin 2X Y 
0 


+ > (n? + 4X") exp (— n?/4)(2X sin 2X Y + msinh nY)}. 


H. E. SALZER 
NBSCL 
This work was sponsored in part by the Office of Air Research, AMC, USAF. 
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8 A. M. TurtnG, “A method for the calculation of the Zeta-function,’’ London Math. 
Soc., Proc., s. 2, v. 48, 1943, p. 180-197. 


®H. G. Dawson, ‘On the numerical value of 4* e7*dx,’’ London Math. Soc., Proc., s. 1, 
v. 29, 1898, p. 519-522. 
10 NBS, Tables of Probability Functions. V. 1, New York, 1941. 
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857[ D].—J. K. Lyncu, “Calculation of (sin x)/x tables,’ Postmaster- 

General’s Department, Research Laboratory Report No. 3238, Melbourne, 

C. 1., (April 1949), C.E. 615/R, Case no. 2266. 7 mimeographed pages 

+ drawing. 

These tables of (sin x)/x were computed to facilitate the calculation of 
the harmonic components of periodic impulses. The function (sin x)/x is 
tabulated for x = 0(1°)180°(5°)1800°; 4D, with first differences. The draw- 
ing contains a graph of (sin x)/x for x ranging from 0° to 1600°, together 
with a straight line converting degrees to radians. The statement in the 
Introduction “‘. . . the interval of tabulations is such that linear interpola- 
tion will suffice,” is not quite correct, since second differences are needed for 
full accuracy when x lies between 180° and 295°. 

For another table of this function, see [MTAC, v. 4, p. 80]. 


H. E. SALZER 
NBSCL 


858[F, G, O].—R. E. Brearp, “On the coefficients in the expansion of e“ 
and e~*.”” Inst. of Actuaries, Jn., v. 76, 1950, p. 152-163. 


The expansions mentioned in the title are written 


ll 


ao 
ec‘ =e > K,t"/n! 


n=0 


ee“ =e") A,t"/n! 
n=0 
where the A’s and K’s are integers. They are tabulated for = 0(1)26. 
[See the following review where K,, is denoted by U(m) and a table is de- 
scribed for n = 1(1)50; the agreement for m < 26 is exact. ] The author used 
the formulas 


K, = > A*0°/k! 
k=1 
and 
A, = > (—1)*A*0"/k! 
k=1 


The values of the Stirling numbers of the second kind, A*0"/k! were taken 
from FISHER & YATES [as corrected in MTAC, v. 4, p. 27] for r < 25. The 
necessary additional values of these numbers for r = 26 are appended on 
p. 163. The table of A, appears to be new. Previous tables of K, for m < 20 
are mentioned in FMR, Index § 4.676. 

D. H. L. 
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s59[F, G, O].—H. Gupta, “Tables of distributions,” East Panjab Univer- 
sity, Research Bull., no. 2, 1950, p. 13-44. 


The author denotes by u(n, a) the number of ways that m unlike objects 
may be distributed into a boxes. If a > m then u(n,a) = 0. The sum 


U(n) = & u(n, a) 
a=1 
is the total number of distributions of m different objects into groups. It also 
represents the number of rhyming schemes in a stanza of m lines. In terms 
of generating functions we have 


= u(n,a)x*/n! = (e —1)*/a!, ¥ Uln)x*/n! = e* exp (€*). 


The fundamental difference equation for u(n, a) is 
u(n + 1,a) = au(n,a) + u(n,a — 1). 


By means of this relation the author has constructed a table (p. 17-43) of 
u(n,a) for 1 < a < nm < 50, and has obtained, by summing, a table of 
U(n) for n = 1(1)50. The function U(n) may also be computed from the 
relation 


(1) U(n + 1) - = uw({): 


This was used to check the whole set of tables by computing the 48-digit 
number U(50). The several terms of (1) for » = 50 are set forth on p. 44. 
D. H. L. 


860[F].—H. Gupta, “‘A table of values of Liouville’s function L(#),”” East 
Panjab University, Research Bulletin, no. 3, 1950, p. 45-63. 


The function of LIOUVILLE, despite the title of the paper under review, is 
usually denoted by \(m) and is +1 or —1 according as the number of prime 
factors of m, each reckoned with its multiplicity, is even or odd. The function 
L(t) is then defined as the sum function 


Li) = & Xn). 
nXt 


According to a conjecture of Potya, the values of L(¢) for ¢ > 1 are never 
positive. This conjecture, when proved, will imply the truth of the famous 
RIEMANN hypothesis. In 1940 the author constructed a table of L(¢) for 
t = 1(1)20000 [MTAC, v. 1, p. 201]. The present paper presents a con- 
densation of this table from which with very little effort an isolated value 
in the original table may be found. More specifically, the present table gives 
— L(t) for t = 5(5)20000 and, by a clever device, the values of \(¢) when ¢ is 
not a multiple of 5. (A(S5k) = — X(R).) 

Two small tables (p. 47) give some idea of the range of values of L(t) 
for t < 20000. The conjecture of Polya is verified but nevertheless —L(¢), 
though positive, is fairly small. A maximum of 150 occurs at ¢ = 15810. The 
function ¢“Z?(¢#) has a maximum at ¢ = 9840. This would indicate that 
L(t) = O(¢*) which would also imply the Riemann hypothesis. The author 
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mentions that isolated values of L(t) up to ¢ = 60000 may be found by 
another method. By still another method the reviewer has found that 
L(4000000) is very close to —1100. 


. mM. 1. 


861[F].—D. JARDEN, “‘On the numbers V5, (” odd) in the sequence associated 
with Fibonacci’s sequence,” Riveon Lemat., v. 4, 1950, p. 38-40 [Hebrew 
with English summary]. 


This note contains a factor table of the numbers 
Von — 5U,4+ 3 and Von + 5U,4+ 3 


whose product, in case is odd, is Vsn/V,. Here V» is the mth term of the 
FIBONACCI sequence: 2, 1, 3, 4, 7, --- 


Ve = Viesd + ¥ nnd IU sm - Vin + aT ask 


The table is for m = 1(2)77. Factorization is incomplete in many cases. For 
the function V2, — 5U, + 3 complete factorization into primes is given only 
for m = 1(2)35, 39(2)45, 49, 51, 63, 75; for the second function only for 
m = 1(2)45. Most of the results come from the table of KRAITCHIK.! 
ee ok. 
1M. Kraitcaik, Recherches sur la Théorie des Nombres. Paris, 1924, v. 1, p. 80. 


862[K].—A. N. Back, ‘Weighted probits and their use,” Biometrika, v. 37, 
1950, p. 158-167. 


The paper tabulates 
Y-5 
P = (2x)-3 f exp (— $u?)du, w= Z(PQ)-, 


wy = w(Y — PZ-), Wym = w(Y + QZ-") 


all to 3D for Y = 1(.02)9, except that wy is replaced by 100 — wyp for 
Y > 6.42. As usual P + Q = 1 and Z = (2x)-* exp (—(Y — 5)?/2). The 
purpose of the table is to minimize the arithmetic of the maximum likelihood 
estimation of the probit regression. There are detailed instructions -for 
calculation. 


H. W. Norton 
University of Illinois 
Urbana, IIl. a 


863[K].—F. N. Davin, ‘‘Note on the application of Fisher’s k-statistics,” 
Biometrika, v. 36, 1949, p. 383-393. 


This paper contains a table of the terms which must be added to «(r*s‘), 
the product cumulant of order hi of the r-th and s-th &-statistics about the 
population 7-th and s-th cumulants respectively in terms of population 
cumulants, in order to give the corresponding product moment. The table 
gives all terms for h + / = 4 and 5 and those for h + / = 6 involving n~* 
only where is the sample size. The paper illustrates the utility of these 
tables in reducing the algebra in finding the approximate moment character- 
istics of certain statistics. 
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864[K].—F. N. Davin, “Two combinatorial tests of whether a sample has 
come from a given population,” Biometrika, v. 37, 1950, p. 97-110. 


In connection with the tests herein described and illustrated, the au- 
N 
thor tabulates (Table 1a, p. 109) N-* € ) A*-20" for N = 3(1)20 and 


Z = 0(1)19 to 4D. This is the probability of Z zero groups when N sample 
units are randomly arranged in N groups of equal probability and is thus 
related to a classical distribution problem. Table 4a, p. 110 gives values of 


N 
(2N)-* fs ) A’-2(N*), for N = 1(1)10 and Z = 0(1)10 to 4D, which is 


the probability that if N balls be dropped at random into 2N compartments, 


Z of a specified set of N compartments will be empty. 
Ss Se On 


865[K].—F. E. Gruss, “Sample criteria for testing outlying observations,” 
Annals Math. Stat., v. 21, 1950, p. 27-58. 


The author proposes the following statistics 











n—1 
s2 2 @-%) 1 In — = 
ee ee ee lee 
a—1 S 
LX (az: — 2)? 
i=1 
hn 2 
se UL & #1) 1 — 2 
3 =1— — Ty, ™% = 
n n—i1 S 
LX (x: — 2)? 
i=1 
h Sa1S ++ S28 = —— Eh = — E28 LE 
where z, < 2 < S fm Bn = Oy to = Bm 


for testing the greatest and lowest observations respectively in a sample of 
size n, from a normal population. Values of S,?/S* or S,°/S* to 4D are given 
at percentage points 1, 2.5, 5 and 10 for m = 3(1)25 in Table I (p. 28) and 
corresponding values to 3D of T, or T; in Table IA (p. 29). Grubbs’ values 
in Table IA calculated from the exact distribution concur with those of 
E. S. Pearson & C. CHANDRA SEKAR! for common entries. 

Table II (p. 31-37) contains values of P(u, < u) to5D for u = 0(.05)4.85, 
n = 2(1)17, u = .5(.05)4.90, m = 18, 19 and to 4D for u = .5(.05)4.90, 
n = 20(1)25, where u, = (x, — #)/o in samples of size m drawn from a 
parent normal population with zero mean and standard deviation ¢. These 
values were obtained from the cumulative distribution 


Fa(u) = 2¥(2e{n — 1})-1 [exp (— © n/(n — 1)) Fa-a(nz/(n — 1))dz 
0 2 


where F,(u) = fo" dF(un) = P(u, < u). This distribution (obtained in 
1945) was also independently arrived at by K. R. Nair? and is equivalent 
to the result of A. T. McKay’ arrived at by means of characteristic func- 
tions. Nair? also gave tables of this probability integral for w = 0(.01)4.7, 











74 RECENT MATHEMATICAL TABLES 


nm = 2(1)9, to 6D which concur in corresponding values with Grubbs’ 
Table II. 

Table III (p. 45) contains values to 3D of u, = (x, — #)/o for percentage 
points 90, 95, 99, and 99.5, m = 2(1)25, also given by Nair? to 2D for per- 
centage points .01, .5, 1, 2.5, 5, 10, 90, 95, 97.5, 99, 99.5, 99.9, m = 3(1)9. 

Table IV (p. 46) contains the mean, standard deviation, a; and a, of u, 
to 4D for m = 2(1)15, and to 3D for m = 20, 60, 100, 200, 500, and 1000. 
These values agree with corresponding values of TrpPpETT‘ for the mean and 
those of McKay’ for the standard deviation except for » = 2, which Grubbs 
gives correctly as 0.4263 (by McKay’s formula) whereas McKay gives 
0.6179. 

Table V contains percentage points 1, 2.5, 5, and 10, m = 4(1)20 to 4D 
for Se-1, n/S? or Si,2/S? where 


n—2 


; a = _ (x; aa 2.20)", Fn-1,2 = (n ya 2) > Li, 


i=1 i=1 


n n 
Siz = X (@i — 412)’, Zi2 = (n —2)") 2, 

i=3 i=3 
S? and & as above. These statistics are proposed for testing the significance 
of the two largest and the two smallest observations in samples of size n 
drawn from parent normal populations. 

S. B. LITTAUER 

Columbia University 
New York 


1E. S. Pearson & C. CHANDRA S:xkar, “The efficiency of statistical tools and a cri- 

terion for the rejection of outlying observations, ” Biometrika, v. 28, 1936, p. 308-320. 
7K. R. Narr, “The distribution of the extreme deviate from the sample mean and its 
studentized form,” Biometrika, v. 35, 1948, p. 118-144. 

*A. T. McKay, “The distribution of the difference between the extreme observation 
and the sample mean in samples of m from a normal universe,”’ Biometrika, v. 27, 1935, 
p. 466-471. 

*L. H. C. Tippett, ‘On the extreme individuals and the range of samples taken from 
a normal population,”’ Biometrika, v. 17, 1925, p. 364-387 


866(K].—A. Hap, ‘(Maximum likelihood estimation of the parameters of 
a normal distribution which is truncated at a known point,” Skandi- 
navisk Aktuarietidskrift, v. 32, 1949, p. 119-134. 


This paper includes four short tables (p. 130-134) of certain functions 
for use in estimating the mean and variance of normal populations from 
one-sided truncated samples with known points of truncation, both for the 
simple truncated case in which the number of missing (unmeasured) observa- 
tions is unknown and for the case which the author labels as “‘censored”’ in 
which the number of unmeasured observations is known. Tables I and II 
pertain to the simple truncated case. Tables III and IV furnish similar 
information for the ‘“‘censored” case. All four tables were computed from 
Fisher’s I or Hh function." 

Table I of z = f(y) is an inversion of Table XVI from the B.A.A:S. 
Tables! where z is the abscissa of truncation in standard units of the com- 
plete distribution and y is a function of moments of the truncated distribu- 
tion. In practice, one sets y = m >> x,2/2(> ~,)? and reads z directly from 
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this table. Entries of z are given to 3D for y = .55(.001).91. In the nota- 
tion of Pearson & LEE,? y = (Wy, —1)/2, and in the notation of FisHER® 
y = I,I2/I:?. Table II contains g(z), a function employed in estimating the 
variance, together with its first differences as an aid in interpolation. Entries 
are tabulated to 4D(4 or 5S). This table also includes elements of the 
variance-covariance matrix, g1:(2), wi2(z), and pe2(z) tabulated to 4S, and 
p(z), the correlation coefficient between estimates of population parameters, 
to 3D. The argument, z, ranges over —3(.1)2. 

Table III of zs = f(h,y) was also computed by Statistika I/S for 
h = .55(.05).8. It permits the abscissa of truncation, 2, for the “censored” 
case to be read directly for given values of h and y, where h is the ratio of 
the number of unmeasured to measured observations, and y is the counter- 
part of the similarly labeled function for the simple truncated case. Entries 
for z are given to 3D with y = .5(.001)1.5 and h = .05(.05).8. The functions 
of Table IV for the “censored” case correspond to those of Table III for 
the simple truncated case. The range and the interval of the argument, z, 
are the same in both tables, and the number of significant digits is the same 
in most instances. 

A. C. COHEN, JR. 

University of Georgia 
Athens, Georgia 

1B.A.A.S., Math. Tables, v. 1, 2nd ed., Cambridge, 1946, Tables XV and XVI. 

? KARL PEARSON & ALICE Ler, tay the generalized probable error in multiple normal 


correlation,’’ Biometrika, v. 6, 1908, p 
3 See 1st ed. of above B.A.A.S. "chine, 1931, p. XXvi-xxxv. 


867[K].—H. O. Hartiey & E. S. Pearson, “Table of the probability inte- 
gral of the ¢-distribution,” Biometrika, v. 37, 1950, p. 168-172. 


This table gives the cumulative probability integral of Student's t for y 
degrees of freedom, namely 


Po,» = {r (24+) /cranvera] f° a + yaya 


to 5D for y = 1(1)20; ¢ = 0(.1)4(.2)8 and y = 20(1)24, 30, 40, 60, 120, ~; 
t = 0(.05)2(.1)4, 5. Also values of ¢ and @ such that 1 — P(t, y) = @ are 
given for a = 10-*, 10-*, 10-5, 5(10)-*; y = 1(1)10. For y = 11(1)20 the 
values of a to 10-5, as well as values of a to 3(10)~* and less for larger y, 
are covered in the main part of the table. Formulas for both single-entry 
and double-entry interpolation are given. 

This table may be particularly useful in problems where it is necessary 
to interpolate for fractional degrees of freedom. It is also useful in cases 
where distributions of other statistics can be expressed in terms of the 
t-distribution. Previous tables' were to 4D for y = 1(1)20. 

F. J. MASSEY 





University of Oregon 
Eugene, Oregon 


1 Student, ‘New tables for testing the significance of observations,” Metron, v. 5, no. 3, 
1925, p. 105-120. 
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868[L].—C. J. Bouwxamp, (I) “On the characteristic values of spheroidal 
wave functions,” Philips Res. Rep., v. 5, 1950, p. 87-90. (II) “On the 
theory of spheroidal wave functions of order zero,’’ Nederl. Akad. 
Wetensch., Proc., v. 53, 1950, p. 931-944. 


I. The differential equation for spheroidal wave functions is taken in 
the form 


(1) (1 — 2*)y” — 2zy’ + (A — m1 — 2°) + R*e*)y = 0. 


Corresponding to a countably infinite set of characteristic values \, for given 
m and k, equation (1) admits a solution which is regular at z = + 1. It is 
known that the form for d is 


Ae (k) = > pik. 
i=0 
The author gives the coefficients p;, 71 = 0(1)4, for general m and m, along 
with numerical values of the coefficients (in both fractional and decimal 
form) for m = 1, m = 1(1)7. 

Three brief tables are given: 

Table 1: Prime factors in the denominators of the coefficients ;, for 
nm = 1(1)7; m = 1. 

Table 2: Decimal form of the coefficients p; form = 1(1)7;m = 1, to 12D. 
MEIXNER (1944) gave a more complete table, to 10D. [MTAC, v. 3, p. 
524-6. ] 

Table 3: Numerical values of \4(k), as calculated from power series, for 
R= +1, +4, n = 1(1)7; to 8D for k? = +1, and 5D for k? = + 4. 

II. In this paper the author treats various series representations for the 
solutions of the first, second, and third kinds, associated with the character- 
istic values of the equation. In particular he derives a set of coefficients c,, 
such that two independent solutions are of the form 


y(s) = exp (— ks) 4 cP.(s);  Y¥*(s) = exp (— is) ry CaQn(2), 


where P,(z) and Q,(z) are the Legendre functions. Here y(z) is the regular 
solution. These solutions are valid in the domain of regularity of the function 
(outside the branch cut joining z = 1 and zg = — 1). Interesting relations are 
obtained between the eigenvalues y,(k) of the integral equation 


y(s) =u f exp (kel) y(t)dt 


and the coefficients cy. Although the series involving c, are important theo- 
retically they do not converge as rapidly as the classical series involving 
Bessel functions. 


GERTRUDE BLANCH 
NBSINA 


869[L].—B. W. Cono tty, ‘‘A short table of the confluent hypergeometric 
function M(a, y, x),”’ Quart. Jn. Mech. Appl. Math., v. 2, 1950, p. 
236-240. 


Tables of confluent hypergeometric functions are listed in section 22.56 
of the FMR Index, and some have also been reviewed in MTAC. In the 
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present paper the author tabulates 


M(a, 7, x) = ry P(y)P(a + r)x"/(r!P (a) (y + 7) = 4 eat 


to 11D for x = .1, .2(.2)1, a = — 1(.2)1, and y = .2(.2)1. The tabulated 
values are believed to be correct to within 2 or 3 units of the last decimal 
place. 


The coefficients a,, and the values of M(a, y, 1) were provided by J. C. 
P. MILLER. 
A. E. 


870, L].—A. Coomss, “The translation of two bodies under the free surface 
of a heavy fluid.’”” Cambridge Phil. Soc., Proc., v. 46, 1950, p. 453-468. 


In the course of the work a 6D table is constructed of 


{li(et) — wijexp (— a(1 — 978)) - Dan! YO (}iaf)*/s! 
n=0 s=n+l1 
for 8 = — 1, 1 and a = 1, 2, 2.5, 3(1)6, 8. More extensive tables are pro- 
vided for the forces acting upon two cylinders in the same vertical or hori- 
zontal plane. 
A. E. 


871[L].—Otto EMERSLEBEN, “Die elektrostatische Gitterenergie eines neu- 
tralen ebenen, insbesondere alternierenden quadratischen Gitters.”’ Zeit. 
fiir Physik, v. 127, 1950, p. 588-609. 


In the study of crystal lattices there appear the EpsTEIN zeta functions, 
defined for R(s) > 2 by 


0 0 


», nl = ee (2ai(kihi + kahz)) 
1 2 


ko hyo (ki? + k;?)*” 
For various integer values of s the author evaluates these functions in the 
four cases (hi, he) = (0,0), (4, $), (0, $), (4, 3). The last three cases can be 


computed directly from the first with the help of functional equations. 
For s = 1 the formula 


0 0 
0 0 


Z 











(1) = —4425 wrn)[1 — F((wn)!)] 


n=1 


Z 








is used, where r2(m) is the number of ways m can be written as a sum of two 
squares and 


F(x) = at f e~"dt. 
0 


Eight terms of the series suffice to yield the value 


0 0 


Zl0 0 








(1) = — 3.900 264 92. 


For s = 2 the values have been previously given by the author.' For s = 3 
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and s = 4 use is made of the product representation 


0 0 
0 0 


and tables of ¢(s), L(s). Here 


Zz (s) = 45(3s)L(4s) 








is) = Sa and = Ls) = ¥ (-1)*"*(2n — 1). 


n=1 n=1 


An alternative method is also given in the case s = 4. For negative integer 
values of s the author uses the functional equation which gives the value at 
s in terms of the value at 2 — s and gamma functions. 


Tom M. Aposto.i 
California Institute of Technology 
Pasadena, California 


1 Otto EMERSLEBEN, “‘Einige Identitaten fiir Epsteinsche Zetafunktionen,”’ Math. Ann., 
v. 121, 1949, p. 103-106. 


872[L].—B. FRIEDMAN, “Theory of underwater explosion bubbles,’ Comm. 
on Pure and Applied Math., v. 3, 1950, p. 177-199. 
Tables are given to 4S of the functions 


flax) = 2x 3 (—1)*[(2n + 1)? — 2} 


F(x) = (1 — x)[$x“f(x) + $n 2 — (1 + x)“]). 
The range for x is 0(.05)1 for f and —1(.05)1 for F. 


873(L].—E. L. Kaptan, “Multiple elliptic integrals,” Jn. Math. Phys., v. 29, 
1950, p. 69-75. 


The author shows that the evaluation of integrals of the form 


f K kk'"'dk, f E k*k'"'dk 


where n, n’ are integers and k, k’, K, E are the standard notations, and the 
evaluation of some similar integrals can be reduced to the computation of 
the eight integrals 


ff Kaj, f Kr ae/e, f Kavye, f Bavx, 
f Ka, ff Kae, f awe, fx dk/k”, 


and that the last four can be dispensed with if Landen’s transformation is 
used in the process of reduction. The values of the 8 basic integrals taken 
between the limits 0 and 1, and modified if necessary to produce convergence, 
can be expressed in terms of known constants and are given in the paper to 
10D. There is, moreover, a 10D table of 


(n + 1) [Kar 
0 
for m = —.9(.1)2. 
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874[L].—A. Leitner & R. D. Spence, “The oblate spheroidal wave func- 
tions.’”’ Franklin Institute, Jn., v. 249, No. 4, 1950, p. 299-321 [see also 
ref. in MTAC, v. 3, 1948, p. 99-101]. 


The differential equations for the “‘angular” and “radial” spheroidal 
wave functions of order m respectively are written in the form 


(1) d{(1 — 9°)du/dn}/dn — {m1 — o°)* — a + (1 — 9) }u = 0, 
(2)  d{(1 + &)dv/dt}/de — {m*(1 + P)* — a+ e(1 + #)}v = 0. 


Here « = ka (k is a given parameter and a is the radius of the focal circle) 
and a is a separation constant whose eigenvalues aj, are such that a solution 
tim belonging to aim is finite at 7 = + 1. 

After a summary of the methods for representing the eigenvalues a;,, and 
the functions #im(n) and vim(£), expansion formulas for spherical and plane 
waves in terms of the spherical wave functions are given. The numerical 
results consist of tables and curves [p. 314-320] and contain for « = 1(1)5 
the following quantities: 


The eigenvalues aim for m=1, 1 = 4(1)8; m = 2, 1 = 2(1)9; m = 3, 
l = 4(1)10; m = 4,7 = 4(1)11. 
The norm Nim for m = 0,1 = 0(1)5; m = 1,1 = 1(1)6. 
1 

The norm Ni» is defined as f Urmtymdn = biNim. 

—1 
dim for m = 0,1 = 0(1)5; m = 1,1 = 1(1)6 
Qim for m = 0,1 = 0(1)5. 
The gim and Qim are constants associated with the radial functions. 
Finally the angular function uim(n) for m = 0,1 = 0(1)5; m = 1,1 = 1(1)6; 
m = 2,1 = 2 and 3, and for values of the argument » = 0(0.1)1. 

F. OBERHETTINGER 


California Institute of Technology 
Pasadena 4, California 


875[L].—BriGITtTE Rapon, “Sviluppi in serie degli integrali ellittici,”” Accad. 
Naz. Lincei, Aiti, Mem., Cl. Sci. Fis. Mat. Nat.,s. 8, v. 2, 1950, p. 69-109. 
The usual notations of elliptic integrals, 


Fle, k) = f "wy, k)dy, Ele, ) = f ” AW, Bay, 


a(e,e,k) = [C+ psint yA, B) Pay 
are used, where 
0< ¢ <r, Oo<F <1, and A(y, k) = (1 — sin? y)!. 


The author gives the expansion of both complete (¢ = $2) and incomplete 
(¢ general) integrals in powers of k?, also in powers of k” = 1 — & (in the 
latter case there are also logarithmic terms). In addition, for the incomplete 
integrals there are expansions in powers of A(¢y, k) — cos ¢, in powers of 
1 — A(¢, k), and for F and E also as a trigonometric series. For the integrals 
of the third kind expansions in powers of k? + p are also provided. 
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In each case the first few coefficients are given explicitly, for the others 
either a general form or a recurrence relation. The region of convergence is 
noted, as is the region of “‘practical convergence” where the series converges 
at least as well as an infinite geometrical progression with ratio }. 

A. E. 


876(L).—Wasao SurpaGakI, Table of the Modified Bessel Functions, Kyushu 
University, Physics Department, April 1946, 132 p. 25.4 X 18.3 cm. 


This table tabulates the modified Bessel functions J,(x), K,(x) or 
Vx In(x), Vx Ka(x). 

The functions Jo(x), I:(x), ---, I2:(x) are tabulated for x = 0(.01)1(.02)5; 
Io2(x) for x = 1(.02)5; the functions x'J,(x) are tabulated for x = 5(.04)25, 
n = 0(1)22. 

The functions K,(x) are tabulated for x = 0(.01)1(.02)5, m = 0(1)23; 
the functions x!K,(x) are tabulated for x = 5(.04)25, m = 0(1)23. 

All functions are given to mostly five significant figures, and to facilitate 
interpolation the first differences, divided by 1004, where h denotes the 
tabular interval in x, are tabulated alongside the functions. The purpose is 
to avoid the computation of » = Ax/h in the interpolation. 

H. E. SALZER 
NBSCL 


877[L].—Toxyo NUMERICAL COMPUTATION BuREAU, Report no. 3. Tables of 

Spherical Wave Functions, 1950. VI. 1. (Keisuh-Renkyukjo) 466, Kyoh- 

doh, Setagaya, Tokyo, Japan. 1950, iv, 77 p. 25.4 X 36.2 cm. 

The text (p. i-iv) of this publication is in Japanese although the title 
page is in English. Various tables occupy p. 2-77. The tables are primarily 
concerned with the Riccati Bessel functions arising in connection with the 
wave equation, as follows: 


F(x) = ($4x)*Tngy(x) = x™**f,(x)/(2m + 1)!! 
Ga(x) = (34x)!J_n4(x) = x-gn(x)(2m — 1)!! 
F,! (x) = (n + 1)x"fa(x)/(2n + 1)! 

G,' (x) = nx-"—'9,(x)(2m — 1)!! 


where (2m + 1)!! = 1-3-5---(2n + 1). 
The first tables are of F,, Gn, Fn’, Gn’ for m = 0(1)6, x = [0(.05)10; 6D], 

5. In limited ranges where F,, F,’ have small significant numbers and G,, 
G,' are too irregular to be interpolated, 6D values of ¢n, Gn, Yn, Vn are also 
given, where 

on = log fa(x) — log (2m + 1)!! 

on = log fa(x) + log [(m + 1)/(2m + 1)!!] 

¥n = log gn(x) + log (2m — 1)!! 

Vn = log gn(x) + log [m(2n — 1)!!]. 


Recurrence formulas may be used for m > 7. Earlier tables of this kind, 
by Doopson and others [MTAC, v. 1, p. 233], were at interval never less 
than .1. 
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The above mentioned tables occupy 53 pages. Then come 6D tables of 
values of Cn, On, Cn, On, Pus Qn, Rn, Sn where 


G, + iF, = C, exp (i{x — 44n + 0,}) 
F,' — iG,’ = C, exp (i{x — 44n + 6,}) 

F, — 1G, = (— Qn + iP.) exp (i{x — $rn}) 
G,' + iF,’ = (R, — iS,) exp (i{x — 32n}). 


These functions are given for y = 1/x in the range 0 to .1 at varying 
intervals .01, .005, .002, and 8?. 
Finally, there are 6D tables of fa, fn, ga, gn for m = 1(1)6, and 
x = [0(.05)m, 6D], &, m varying from 2 to 5.35. 
R. C. ARCHIBALD 
Brown University 
Providence, R. I. 


MATHEMATICAL TABLES—ERRATA 


In this issue references have been made to Errata in RMT 865 (Grusss). 


181.—R. L. ANDERSON & E. E. Houseman, Tables of Orthogonal Polynomial 
Values Extended to N = 104 [MTAC, v. 1, p. 148-150]. 


A recalculation of this table reveals the following complete list of errata 


Page Line n For Read 
615 broken type —585 
618 496388 4496388 
— 3583 — 2583 
625 +32481 
629 +2648 
642 135050 
649 +701935 
653 +88686 
663 —2107 
665 +110308 
669 — 26592 


The last of these was reported by W. F. Brown Jr. in MTAC, v. 4, p. 222. 
The values given for » = 72 on p. 640, line 36 are really for » = 71. 
The correct values are: 124392, 966 52584, 15878 63880, 39 89066 92520, 
3436 29622 27080. 

JAcK SHERMAN 


Beacon Laboratories 
Beacon, N. Y. 


182.—N. G. W. H. BEEGER, “On the congruence (p — 1)! = —1 (mod #*),” 
Messenger Math., v. 51, p. 149-150, 1922. 


On p. 150, p = 239 for w = 147 read 107. 
This error was discovered as the result of recalculation and extension 
of Beeger’s table of Wilson’s quotient to p < 1000 by use of the SEAC. 
K. GOLDBERG 
NBSCL 
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183.—H. T. Davis, Tables of the Higher Mathematical Funcztions. V. 2, 
Bloomington, 1935. 


Table 38 of log E,, p. 298 


n for read 
11 6358 6400 
37 9184 9084 
42 0201 0301 
44 3908 3907 


R. LIENARD 
32 Boulevard Saint Michel 
Paris, France 


184.—R. E. GREENWoop & J. J. MILLER, “Zeros of the Hermite poly- 
nomials and weights for Gauss’ mechanical quadrature formula,’”’ Amer. 
Math. Soc., Bull., v. 54, 1948, p. 765-769 [MTAC, v. 3, p. 416]. 


The above table was compared with our similar 13-place manuscript 
table. Apart from a few discrepancies of only a unit in the last place, the 
following errors were found in the weight factors: 


page n for read 
768 7 As = 0.000 971 781 258 Az = 0.000 971 781 245 
10 As = 0.001 343 645 77 As = 0.001 343 645 75 
H. E. SALzer 
R. ZUCKER 
R. E. Capuano 
NBSCL 


185.—Z. Kopat, ‘‘A table of the coefficients of the Hermite quadrature 
formula,”’ Jn. Math. Phys., v. 27, 1948, p. 259-261 [MTAC, v. 3, p. 473]. 
On page 260, corresponding to m = 18, the sixth value of p, 

for 0.000051 59 read 0.000051 80 
H. E. SALzeErR 
R. ZUCKER 
R. E. Capuano 
NBSCL 


186.—M. Kraitcuik, “On the divisibility of factorials,’’ Scripta Math., v. 
14, 1948, p. 24-26 [MTAC, v. 3, p. 357-8]. 


P. 25, m= 18, for 108514808571661 read 226663-478749547 


[The fact that this large “‘prime”’ factor of 18! —1 is composite was proved 
bv A. Ferrier, letter of 16 Nov. 1950, who showed that (for V = 108- - -661) 


2%—! = 5107 05663 62894 (mod JN). 
The problem of factoring N was proposed during a demonstration of the 


SWAC on 29 Jan. 1951. The factor 226663 was discovered by the SWAC 
in 25 minutes running time. D. H. L.] 
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187.—NBSMTP Table of Natural Logarithms. V. 1, New York, 1941. 


P. 184, x = 18254 for 9.18213 read 9.81213 


LINDLEY M. WILson 
NBS 


188.—WILHELM MaGnus & Fritz OBERHETTINGER, Anwendung der ellip- 
tischen Funktionen in Physik und Technik. Berlin, 1949 [MTAC, v. 4, 
p. 23]. 
P. 113, col. 4, last line. For 1,26928 read 1,29628 
T. J. Hiceins 
D. K. REITAN 


Univ. of Wisconsin 
Madison 


189.—A. Reiz, “On the numerical solution of certain types of integral 
equations,” Arkiv Mat., Astr. Fys., v. 29A, no. 29, 1943, 21 p. [MTAC, 
v. 3, p. 26.] 


On p. 6, Reiz tabulates +x;, p; and a;, where x; are zeros of the Hermite 
polynomials, p; are x~-weight factors, and a; are weight factors-exp(x,), 
for m = 2(1)9. The following errors of more than a unit in the last (7th) 
decimal place were found: 





for read 
2 = 3, 1.3239316 1.3239312 
n= 4, 1.2402244 1.2402258 
n= 5, 1.1814877 1.1814886 
n = 6, 1.3358489 1.335849] 
0.9355808 0.9355806 
1.1368912 1.1369083 
2 = 7, 0.8971839 0.8971846 
1.1013979 1.1013307 
n = 8, 1.9816821 1.9816568 
0.866838 1 0.8667526 
1.0718011 1.0719301 
n= 9, 0.8417403 0.8417527 
1.0449691 1.0470036 
H. E. SALzer 
R. ZUCKER 


R. E. Carpuano 
NBSCL 
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[Eprtor1aL Note. The Mathematical Tables Committee of the Royal 
Society wishes to announce the establishment of a depository of unpublished 
mathematical tables in the library of the Royal Society. Lists will be pub- 
lished periodically of the tables which have been accepted. The tables will 
be available for examination in the library and it is proposed to arrange for 
photo-copies to be supplied as a reasonable charge to those who desire them. 
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Communications should be sent to the Assistant Secretary, The Royal 
Society, Burlington House, Piccadilly, London W.1]. 


116[C, F].—NBSCL—Radix table for the computation of logarithms to many 
places using logarithms of small primes. Preliminary manuscript awaiting 
final check in possession of NBSCL and deposited in UMT Fite, 45 
leaves, hectographed. 


This table, which was calculated at the suggestion of J. B. Rosser, gives 
the integers a, 5, c, d, e, f, and g and the corresponding quantity 


N = 293°5°7411°13/179 


for about 2000 values of NV ranging from 9.996 to about 1 + 2-10-"*. Exact 
values of NV are given down to NV = 1.155, thereafter N is given in the form 
1 + v, where v is given to 5S. 

The purpose of this table is to facilitate the computation of logarithms to 
very many places, making use of the fact that the logarithms of 2, 3, 5, 7, 11 
and 17 have been given by H. S. UHLER! to 330 decimals. The logarithm of 13 
was calculated to 328 decimals, using Uhler’s values. The logarithm of 37 is 
obtained as a by-product. 

The method of construction and the method of using the table are de- 
scribed in the introduction. The entire computation and most of the pre- 
liminary checking was done by ELIZABETH GODEFROY, under the super- 
vision of H. E. SALZER. 

1H. S. Unter, ‘Recalculation and extension of the modulus and of the logarithms of 


2, 3, 5, 7, and 17.” Nat. Acad. Sci. Washington, Proc., v. 26, 1940, p. 205-212; “Natural 
logarithms of small prime numbers,” Ibid., v. 29, 1943, p. 319-325. 


117[F].—BALLIisTIC RESEARCH LABORATORY, Table of Fermat's quotients, 236 
leaves, tabulated from punched cards. On deposit in UMT Fine. 


Let p = p, be the n-th prime and let e = e(p,) be the least positive inte- 
ger x such that p divides either 27 + 1 or 2* — 1 (in other words ¢(p) is the 
exponent of 4 modulo ). Finally let r() and R(p) denote the least non- 
negative remainders on division of 2‘ — 1 by p and #? respectively. Then 
the main table gives 


n,p= Pn; P’, “= e(p), f it (p os 1)/e, r(p) and R(p) 


for n = 2(1) 2860, that is for p < 26000. This table was computed in 1949 
on the ENIAC by GreorGE REITWIESNER and Homé MCALLISTER during 
check periods at the suggestion of D. H. L. and is intended to give informa- 
tion about Fermat’s quotient! 


gz = go(p) = (2?-* — 1)/p. 
This integer is related to R and f by the congruence 
gz = $fR(R + 2)/p (mod ). 


The interesting case of g2 being divisible by is thus equivalent to R being 
0 or p — 2. This occurs only twice in the whole table, at p = 1093 and 
p = 3511. These exceptional primes were discovered by MEISSNER (1913) 
and BEEGER (1922); Beeger’s calculations? are for p < 16000. The main 
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table has been rearranged in three other ways, by sorting the cards on which 
the results were punched, as follows 


1) In two parts according as 2‘ = 1 or —1 (mod p) 

2) According to values of f 

3) According to values of f in each of the cases 2‘ = 1 or —1 (mod ). 

1 See L. E. Dickson, History of the theory of numbers. V. 1, Washington, 1928, Chapter 4. 

2 N. G. W. H. BEEGER, ‘“‘On a new case of the congruence 2?-' = 1 (mod p*),”” Messenger 
Math., v. 51, 1922, p. 149-150. 


“On the congruence 2?! = 1 (mod #*) and Fermat's last theorem,” Nieuw Archif v. 
Wiskunde, s. 2, v. 20, 1939, p. 51-54. 


118[F].—A. GLopEN, Tables des solutions des congruences X*" + 1 = 0 (mod p) 
n = 4, 5, 6; p < 10‘. Typewritten manuscript, 8 leaves. Deposited in 
UMT Fite. 


All solutions <p/2 of the congruences 
x = —1, 2! = —1, x*4 = —1 (mod p) 


are given for those primes p < 10‘ for which such congruences are possible. 


119[F].—A. GLopEN, Tables de décomposition des nombres premiers 8k + 1 
dans l’intervalle 350000—600000 en a? + 6* et 2c? + d’ et tables des solutions 
des congruences ? +12=0 et 2m?+12=0 pour ces mémes nombres 
premiers. Handwritten manuscript of 97 leaves, deposited in UMT Fie. 
The original manuscript deposited in the Bibliothéque National de 
Luxembourg. 


This is an extension of the author’s UMT 113 [MTAC, v. 5, p. 28]. The 
author plans a further extension to 800000. 


120[F].—R. M. Rostnson, Stencils for the solution of systems of linear con- 
gruences modulo 2. Box of 1024 IBM punched cards on deposit in the 
UMT Fie. Copies of these stencils may be obtained at cost from the 
Computation Laboratory, University of California, Berkeley 4, California. 
Any system of 9 or fewer linear congruences may be solved modulo 2 by 

simply superposing and sighting 9 or fewer punched cards selected from the 

set. This set of stencils is used in the application of a sieve method for the 
solution in integers of systems of linear equations devised by H. Lewy. 


121[K].—D. L. GrcBert, Levels of significance—a direct table of the two-tail 
probabilities. 38 leaves, mimeograph manuscript deposited in the UMT 
FILE. 

This is a direct tabulation of the areas in the tails of the Student-Fisher 
distribution, that is, of 


t 
P= P(|u|>t) = 1 — n-*(B(jn, yf (1 + u2/n)-@Pdy, 
-t 
Values of 100P are given to 1D for m = 1(1)35 and ¢ = 0(.01)5.99. There is 


a supplementary table for m = 1, 2, 3 which gives the ranges of values of ¢ 
to 2D for which P(u = #) = 0(.1)5.5. Previously published tables of direct 
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values (most tables, including the recent ones, give values of ¢ for which 
P(|u| = #) has selected values) have proceeded by steps of .1 for ¢, the two 
principal ones being due to W. S. Gossett (‘‘Student’’) [ Metron, v. 5, No. 3, 
1925, p. 105-108], which in its main part gives values to 4D for m = 1(1)20 
and ¢ = 0(.1)6, and to Kart Pearson [Tables for statisticians and bio- 
metricians, Part 1, 3rd edition, London, 1930, p. 36]. The present tables 
were calculated by linear interpolation for 0 = ¢ = 6 and 1 = » = 20 from 
Gossett’s 1925 tables; for higher values of ¢ and m, approximate methods 
given by Gossett were used. [Some spot checking by direct calculation of 
the corresponding incomplete 8-functions for ¢ = 1(.01)1.09 with m = 1 and 
10 gave in 8 cases values differing from Gilbert’s by 1 in the third decimal 
place. For a more critical value » = 20, ¢ = 0.5 a direct evaluation of the 
incomplete 8 integral confirmed Gilbert’s value of 0.623. C. C. C.] 
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Cannon, 225 Far West Building, National Bureau of Standards, Washington 25, D. C. 


TECHNICAL DEVELOPMENTS 


Our contribution under this heading, appearing earlier in this issue, is 
“The California Digital Computer,” by Pau L. Morton. 


DISCUSSIONS 


“Floating Decimal” Calculation on the IBM Card 
Programmed Electronic Calculator 


There is a wide variety of problems occurring in run-of-the-mill com- 
puting in which it is extremely helpful to represent numbers x as xo-10? 
where x» has the same significant digits as x within the capacity of the 


machine used, 1 < xo < 10, and # is an integer. Three types we might 
mention are: 


1) Problems in which certain computed quantities have magnitudes 
difficult to estimate. 

2) Problems in which one or more quantities have such a wide range of 
magnitude that no single fixing of the decimal point will suffice for 
the entire range. 

3) Problems in which the setup time involved in estimating magnitudes 


is not justified by the saving in machine time made possible by fixed 
decimal calculation. 


For application of the IBM Card Programmed Electronic Calculator to such 
problems, we have devised a calculator programming of the general-purpose 
type based on the “floating decimal” representation of numbers. A brief 
explanation of the way the machine is instructed will help in an understand- 
ing of what follows. 

The calculating unit (IBM 604) of the Card Programmed Electronic 
Calculator is, for one step in the calculation, instructed to receive numbers 
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A and B, to perform on them an operation f;(A, B), and to leave result C 
where it can be used for further calculation or can be stored. The numbers 
A, B, and C may have up to ten digits and a sign. The choice of operations 
f; is determined to suit the needs of the group using the machine. The more 
calculation one can put into one operation, the faster the C.P.E.C. can do 
the problem. Hence, for a large problem, it is sensible to design the operations 
f; with this in mind. A more versatile but slower setup is one in which the 
same set of operations will suffice for any problem. It will rather obviously 
include among its operations addition, subtraction, multiplication, and 
division. What additional operations are included is up to the needs of the 
group using the machine. The choice of operations is put into effect by wiring 
of plugboards. Most operations are completely instructed by wiring of the 
plugboards for the 604, but on occasion, it is useful to alter A or B before 
they are delivered to the calculator or C as it comes from the calculator. This 
implies special wiring on the tabulator plugboard. 

In our programming of the C.P.E.C. for general-purpose floating decimal 
calculation, we represent a number x by a ten-digit number X. In the repre- 
sentation x = x9-10?, we carry eight significant digits in x» and allow p the 
range — 50 < p< 49. The first eight digits of X are the digits of x». The 
last two digits of X are those of p + 50. The sign of X is that of x». Ex- 
amples are: 


x xX 
4.1752139- 104 4175213954 
— 1.9920367 - 10-*° — 1992036730 
The set of operations which is available in our setup is as follows: 


fo(A,B) =A 
fi(A,B) =A+B 
f(A,B)=A—B 


f(A, B) = A-B 
f(A, B) = A/B 
f(A, A) = VA 


f(A, B) =|A|-B 
fA, B) =|A\|/B 


Corresponding toi = 3,4, 5, 6, 7 there is an operation f’; such that f’; = —f;. 

If no number is delivered to the calculator to be used as A or B ina 
calculation, then the number already in that position is used. Because num- 
bers delivered to storage cannot be read from storage for use as A or B 
immediately, it is often convenient to store the result of a calculation in the 
A or B position in the calculator before the next instructions are given. 
For all operations except fo, we automatically store C in the A and B posi- 
tions. For fo, C is stored in the B position, and the unused B in the A position. 

Operations 6 and 7 are made possible by subsidiary wiring on the tabu- 
lator plugboard which replaces A by —|A]|. Hence, neither of these opera- 
tions can make use of holding the previous result within the calculator in 
A position. 

Our floating decimal programming for the calculator relies heavily on 
the program repeat feature of the 604 used in the C.P.E.C. This feature 
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allows the machine, given a set of instructions, to repeat it several times in 
one operation and, if desired, to use one subset of these instructions on one 
time through a set, a different subset the next time. This use of program 
repeat will make addition or subtraction take two or three times the normal 
computing time if A and B are different by several orders of magnitude or 
if C has a magnitude much smaller than that of A and B. As is to be expected, 
the square root operation will always be slow. The other operations will 
always be finished in the minimum time for a single operation. 

Although it is rather expensive in selectors on the tabulator plugboard, 
one can easily: instruct the machine to make a reasonable (within a factor 
of 4) guess at the gth root of a number which then is improved by iteration. 


: [? +50] —50 
The guess at the gth root of xo- 10? is 3-10‘? 


integer <x. 

The zero of our “floating decimal’’ setup cannot conveniently have all 
its digits zero. An explanation of our addition will demonstrate the origin 
of our representation. Suppose two numbers are added which differ in abso- 
lute value only in the eighth significant digit but which are of opposite sign. 
Then the digits of the result will lead off with zeros. In order to get the answer 
back to standard representation, we examine the left-hand digit. If this 
digit is zero, we shift the digits of the answer one position to the left and 
decrease the calculated p by one. This shifting is continued until the left- 
hand digit is non-zero. To enable the machine to stop shifting if all eight 
digits of the answer are zero, we subtract one in the right-hand digit of the 
answer each time a shift is made after the first. Thus on our machine, 


where [x] is the greatest 


xo°10? — xo-10® = —1.1111111-107. 


This shifting device gives an obvious means for converting fixed decimal 
numbers of eight digits or less to floating decimal numbers. Suppose the 
number of largest magnitude in a column of figures is x-10?. We punch 
x10? into the card with standard “floating decimal” representation and all 
the other numbers with the same last two digits and the same decimal posi- 
tion. Then transferring these numbers through the 604 will convert them 
to standard “floating decimal” representation. 

One additional operation which is useful only in the printing or punching 
of final results has been included in our planning. If calculations made on 
the C.P.E.C. are to be used on other IBM machines, it may be desirable to 
convert results to the fixed decimal form. This procedure will often replace 
key punching by reproducing. The board is so arranged that if A = x94-10"4, 
B = xo®-10"2, pp > pa, and the operations code is 2 and 7, then A will come 
out on Channel C with exponent pz and the left-hand ps — pa digits will 
be zeros. The digits of x9? do not affect this rounding of x94; thus it is not 
necessary to save an entire storage unit to carry pp. 

When Channels A and B are transmitted to the 604, the significant digits 
of Channel A are entered in General Storage 1 and 2, its exponent in Factor 
Storage 4; the significant digits of Channel B are entered in Factor Storage 
1 and 2, and its exponent in Factor Storage 3. The sign entries of FS 1 and 2 
and GS 1 and 2 are wired from the channel sign hubs, in normal fashion. 
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The sign wiring of FS 3 and FS 4 is unusual and important to the opera- 
tion of the board. Its effect is that these two units always read true figures 
with negative sign. Two pilot selectors are used. The first is picked up when- 
ever a spread read-in card is at the third brushes. The second is picked up 
immediately by the impulse from ‘Transfer and S.P. X Control Plus” hub 
for counter 2A, which is wired to total on all cycles. FS 3 sign hub is wired 
to common of one position of the first selector. The transferred hub of that 
selector is wired to the 10 impulse emitter, and the normal hub is wired to 
common of one position of the second selector. Transferred of the second 
selector is wired to the 10 impulse emitter, normal is wired to the sign of 
Channel B. FS 4 is wired similarly in the other position of each of the two 
pilot selectors, with the final connection to sign of Channel A. 

The significant digits of Channel C are in GS 3 and 4, and the exponent 
in the two lowest order positions of the 604 counter. Channel C sign is to 
be taken from GS 3 and 4. 

The absolute value device is this: an x-punch is used to instruct the 
machine; it picks up a pilot selector, which transfers the sign of GS 1 and 2 
from Channel A sign hub and connects it to the sign of FS 4. Then GS 1 
and 2 contains the negative of the absolute value of A. 

The coding that has been used for this board is: 


No Instruction Addition (fi) 
2 Subtraction (f2) 
3 Multiplication (fs) 
4 Division (fs) 
4and 5 Square Root (fs) 
7 Transfer (fo). 


The effect of instruction 2 is just to substitute —B for B. Therefore we also 
have: 


2 and 3 Negative Multiplication (f’;) 
2 and 4 Negative Division (f’s) 
2, 4, and 5 Negative Square Root (f’s). 


If we use the x-punch for absolute value in addition to operation codes, 
we have: 


x, 2, and 3 (fe) 
x, 2, and 4 (fz) 
x and 3 (f’s) 
x and 4 (f’7). 


Calculate selectors 2, 3, 4, 5, and 7 may be picked up in the ordinary 
way. Selector 1 is to be picked up whenever either operation code 3 or 4 is 
present. Selector 6 is to be coupled to selector 5. Selector 8 is to be coupled 
to selector 4. 

On the programming chart, some positions are marked with an asterisk. 
These are positions where the instruction must be selected according to 
which operation is being performed. The selection is indicated in the “‘re- 
marks” column. In that column the symbol 012, for instance, refers to 
program 01, exit 2. 
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The significance of the symbols at the program positions is determined 
by the heading of the column in which they appear. For instance, GS 1 and 2 
in the Exit 1 column means General Storage 1 and 2 read-out. The same 
symbol in the Exit 2 column means General Storage 1 and 2 read-in. The 
abbreviations R & R and RO refer to the counter. We have used the sym- 
bol # to help distinguish between abbreviations for General Storage and 
those for Group Suppression. Thus we have GS 3 and 4 (Generali Storage 
3 and 4) and GS #4 (Group Suppression number 4). We have written mis- 
cellaneous orders (Balance Test for Step Suppression, group suppression 
controls, etc.) in Exit 3 where possible. 

DonaLp B. MACMILLAN 


RIcHARD H. STARK 
Box 1663 
Los Alamos, New Mexico 


APPENDICES 


Suppression Types for General-Purpose Floating Decimal 
Programming for C.P.E.C. 604 


Types of 
suppression: 
I Group suppress # 1 
II Grotp suppress #2 
III Group suppress #3 
IV Group suppress #4 
Vv Suppress on negative balance. 
VI Suppress on positive balance. 
VII Suppress on non-zero balance. 
x Suppress if selector 5 is normal. 
XI Suppress if selector 5 is up. 
XII Suppress on negative balance; suppress if selector 5 is normal. 
XIII Suppress on negative balance; suppress if selector 1 is up. 
XIV Suppress on negative balance if selector 1 is normal; group suppress #1 
if selector 1 is up. 
XV Group suppress #1 if selector 4 is up. 
XVI Suppress on negative balance if selectors 1 and 5 are normal; suppress 
on positive balance if selector 5 is up; group suppress #1 if selector 1 
is up and 5 is down. 
XVII Suppress if selector 4 is up. 
XVIII Group suppress #1; suppress if selector 5 is up. 
XIX Group suppress #1; suppress if selector 4 is up. 
XX Group suppress #4 if selector 5 is normal. 
XXI Group suppress #1; suppress if selector 4 is normal. 
XXII Suppress on zero balance if selectors 3 and 7 are normal; suppress if 
selector 3 is up and 7 is normal; group suppress #1 if selector 7jis up. 
XXIII Suppress on zero balance; suppress if selector 1 is up. 
XXIV Suppress on negative balance; suppress if selector 7 is up. 
XXV Suppress if selector 5 is up; group suppress #2 if both selector 2_and 
selector 7 are up. 
PROGRAMMING CHART 
Read Out Read In Shift Type 
Exit 1 Exit 2 Exit 3 Sup. Remarks 
1 FSi1&2 CTR-+* GS #4 PU II 012 CTR —_ if sel 2 is up 
2 RO MQ XVII 
3 R&R FS1&2 " II 033 out of 6 if sel 3 is up 
4 FS3 CTR -* XVIII 042 CTR+ = if sel 3 is up 
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ned PROGRAMMING CHART—Continued 
id 2 Read Out Read In Shift Type 
me Exit 1 Exit 2 Exit 3 Sup. Remarks 
The 5 FS4* CTR +* ° none 051 FS3 if sel 7 is up and 
7m- sel 2 is down 
and 6 Zero test XXIII 052 «+ if sel 5 is up 
| 7 R&R® CIR-* into? XI 053 B.T.SS. if sel 7 is down or 
age sel 2 is down 
nis- 8 FS4 CTR + XIII O7iemit5 if sel1isup 
sion 9 FS3 FS 4 XIII 072CTR+ if sel 3is up 
10 R&R FS 3 XIV 
ii FS3 CTR + B.T.SS. x 122x+ if sel 3 is up 
12 GS1&2 CTR+* 5 XVI = 123 into 5 if sel 4 is up 
13 GSi1&2 CTR+ into 4 XII 142MQ if sel 3 is up 
14 FS1&2 GS1&2* XIV 142 + if sel 4 is up 
15 R&R FS 1 & 2* ° XIV 152GS3&4 if sel 4is up 
16 FS1&2 CTR+ XXIII 153 out of 6 if sel 3 is up 
17 Half adj. XXIII 
18 R&R FS 1&2 out of 2 XXIII 
19 FS3 CTR + XXIII 
20 Emit 1 CTR — XXIII 
21 RO FS 3 XXIII 
22 FS4 CTR — XXIII 
23 Zero test XXIII 
24 R&R* CTR + into 6 XV 241 GS 3&4 if sel 4is up 
25 MQ GS 3&4 into 2 XXI 262x+ if sel 3 is up 
26 GS1&2 CTR-+* into 3* XIX 263 into 1 if sel 1 is up 
27 FS1&2 CTR-+* into 3* I 272 + if sel 4 is up 
28 GS1&2* CTR — into 3* XXII 273 into 1 if sel 1 is up 
29 Emit 1 CTR — into 3 XxX 
30 GS3&4 CTR+ into 6 XXI 281 R&R _ ifsel4isup 
31 MQ CTR + into 2 XXI = 283 into 1 if sel 1 is up 
32 Halfadj. GS#3PU _  into2 I 
33 R&R GS 3&4 out of 3 x 
#1 34 Emits MQ into 3 none 
35 FS1&2 x+ GS #3 DO xX 
36 GS3&4 - CTR +* into 4* Ill 362 x + if sel 5 is up 
ees 37 RO GS3&4  outof3 XI 363 into 1 if sel 5 is up 
‘or 1 38 Half adj. into 3 none 
39 R&R FS1&2 out of 4 none 
40 Emit 5 FS 3 into 2 x 
41 Emiti* CTR+* GS #4 DO Ill 411 FS4 if sel 5 is up 
42 Emiti CTR — GS#2DO XIX 412x*+ if sel 5 is up 
43 FS3 CTR + XI 
: 44 FS3 s- x 
ss if 45 R&R FS 3 ° none 453 outof4_ if sel 5 is up 
yup. 46 FS1&2 CTR-—- GS #1 PU VII 
47 GS3&4 CTR+ GS #1 DO x 
: 48 RO GS 3&4 out of 3 x 
and 49 R&R B.T.SS. none 
50 MQ CTR — into 5 XXV 
51 MQ CTR — into 5 XXV 
52 FS1&2* CTR+ VI 521 GS 3&4 if sel Sis up 
53 FS1&2* CTR— Vv 531 GS 3&4 if sel 5 is up 
54 
55 R&R GS #3 DO B.T.SS. none 
56 GS #2 PU Prog. Rpt. VI 
57 FS1&2 GS1&2 XXIV 
58 FS3 FS 4 XXIV 
59 FS3 CTR + V 
60 FS1&2 GS3&4 V 
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DIAGRAM OF THE CALCULATOR SELECTORS 
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Notes on Numerical Analysis—4 
Relaxation and Step-by-Step Methods 


Classical methods for the numerical solution of ordinary differential 
equations are all of the step-by-step type, directly applicable only to equa- 
tions with conditions specified at one end of the range of integration. The 
comparatively recent development of relaxation methods provides a new 
technique for the solution of differential equations for which conditions are 
given at two distinct points in the range of integration. We shall describe 
a combination of the two methods which has proved effective in certain 
cases, for instance, when the behavior of the solution is exponential in one 
part of the range and oscillatory in another. 

To fix ideas consider a second-order equation, from which the first 
derivative has been removed, of the form 


(1) y’ +fy = 2, 


where f and g are functions only of the independent variable x. We can 
establish a finite-difference equation to replace (1), connecting values of y 
at three adjacent pivotal points, separated by an interval h, of the form 


(2) rot + Yr—-1 — (2 — Wf.) 9, — hg + A(yr) = 0. 
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Here, A, the difference-correction, is given by 
(3) A = — 6/12 + 66/90 + ---. 


(We could use a more powerful three term formula, involving a much smaller 
A, but these equations are adequate for our present illustrative purpose.) 

An equation of type (2) is to be satisfied at every pivotal point in the 
range of integration. These equations, moreover, are quite independent of 
the initial or boundary conditions of the problem, but the latter will deter- 
mine the method by which such equations are solved. 

If the conditions are of the boundary-value type, we shall commonly 
have specified values yo and y, at the two ends of the range of integration. 
With initial neglect of A, we can therefore solve equations (2) as a set of 
linear simultaneous equations. 

If the conditions are of the initial-value type, we shall be given, or can 
find fairly easily, two adjacent values yo and y, at the start of the range. 
Again, neglecting A, we can in this case solve equations (2) by recurrence, 
obtaining in succession values of ye, ys, ---, etc. 

In each case the approximate solution so obtained can be differenced, 
an approximation to A calculated at every point, and with these values 
inserted as additional terms in equations (2), the respective processes can 
be repeated to obtain better approximate solutions y. This cycle is carried 
on until the equations (2) and (3) are accurately satisfied. 

Now the linear equations are most easily solved when the coefficients 
are dominant in the principal diagonal. The use of iteration or relaxation 
methods is then particularly convenient, and the results are relatively free 
from error. This dominant diagonal is present whenever the solution is of 
the exponential type, when f is negative. If the solutions are oscillatory, 
with f positive, the linear equations are relatively ill-conditioned and their 
solution much more difficult and less precise. 

The recurrence method has precisely the opposite features. In exponen- 
tial regions there is an uncontrollable growth of error; in oscillatory regions 
the error is comparatively negligible. 

We would therefore like to use relaxation in exponential regions, recur- 
rence in oscillatory regions. In some problems of practical interest the solu- 
tion has different forms in different regions, and sometimes the given condi- 
tions allow us to use the appropriate methods. 

Consider, for example, an equation of the form 


(4) y’+Al(x)y=0, (0) = 0, 


in which A(x) is large and negative near the origin, later increasing, passing 
through zero, and having an asymptotic value of say unity. The asymptotic 
solution, for large x, then has the form 


y ~ asin x + bcosx, 


and the problem is to determine the ratio a/b. 
The solution begins by increasing rapidly, then decreases and begins to 
oscillate, as shown in figure 1. 


The finite-difference expression for (4) is given by 
(5) Vr+1 + Se (2 _ h®A,)y, + A(y,) = 0. 
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Now the coefficient of y, is initially greater than 2 in absolute value, and 
begins to be less than 2 when A(x) passes through its zero. Let the point at 
which this occurs be called x,. The solution is exponential for x < x-, oscilla- 
tory for x > x. 

Equation (4) is homogeneous in y, so that y can be given any arbitrary 
value at one particular point other than the origin. Fix the value y = y, 
at x = x,.. For x < x, we then have a boundary-value problem easily soluble 


1 ee TP oS 


ie) Xe — a a 





Fig. 1. 


by relaxation methods. Let us by this means obtain the first approximation, 
with A neglected, in this range. This solution provides at least two starting 
values, at x, and x,.,, from which we can carry forward the solution by 
recurrence into the oscillatory region. 

Since the same equation (5) has been used throughout, the approximation 
so obtained is ‘smooth’ throughout the entire range. We can therefore 
difference and calculate A everywhere. A better solution can now be obtained 
in the exponential region by relaxation. The improved value at x._,, together 
with the given value at x, and the difference-corrections, enable us again to 
carry this better approximation into the oscillatory region, as far as neces- 
sary. The process can be repeated until equations (5) are everywhere 
satisfied. 

There is here a nice computational point. The exponential region is com- 
pletely self-contained, and the final solution, with A included, could be 
obtained before ever proceeding into the oscillatory region. If this were done, 
however, the accurate exponential solution would not fit well with the first 
approximate oscillatory solution: there would be some sort of discontinuity 
at x,, the differences would not be smooth, and the calculation of A uncertain. 
For this reason the application of A is deliberately deferred until it can be 
effected smoothly over the whole range. 

The same technique will help us to throw light on a matter hitherto 
obscure in relaxation. From a first approximation to the solution of a 
boundary-value problem the difference correction is not immediately cal- 
culable near the ends of the range, since central differences do not exist there. 

The first approximation, however, enables us to extend this approximate 
solution, by recurrence methods, outside the chosen range in both directions 
(provided, of course, that the function exists outside this range). We can 
then compute A not only inside, but also at points outside the range. The 
internal values enable us to improve our internal solution, and these new 
internal values, together with the external A’s, allow us to obtain a new solu- 
tion externally, to the same degree of approximation as the new internal 
values. The process is repeated until there is no further correction. The 
number of values determined outside the range depends, of course, on the 
order of significant differences in A. 
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If the solution is everywhere of exponential type, the external recurrence 
may not be very accurate, the error at a point growing with its distance from 
the particular end point. Since these values are required, however, only for 
the computation of A, which has small coefficients multiplying the differ- 
ences, some error can be tolerated. 

L. Fox 
NBSCL and 
National Physical Laboratory 
Teddington, England 
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A brief description of General Electric’s OMIBAC (Ordinal Memory 
Inspecting Binary Automatic Calculator), a binary machine having a 
magnetic cylinder memory particularly useful for repetitive problems. 


2. E. C. BERKELEY, “Simple Simon,” Scientific American, v. 183, Nov. 1950, 
p. 40-43, illustrs. 


The author describes a simple digital relay computer called Simon which 
was constructed to illustrate the principles of automatic digital computers 
[MTAC, v. 4, p. 239]. 

M. M. ANDREW 
NBSMDL 


3. Epmunp C. BERKELEY, ‘‘The relations between symbolic logic and large- 
scale calculating machines,”’ Science, v. 112, p. 395-399. 


This is an elementary discussion of some of the logical problems arising 
in constructing and programming for a modern automatic computer. A 
representative machine, capable of certain simple logical and arithmetical 
operations is described, and the “program” for calculating square roots on 
it is discussed in some detail. 
MICHAEL MONTALBANO 
NBSMDL 


4. J. L. Hix, Design features of a magnetic drum information storage system, 
Engineering Research Associates, Inc., St. Paul 4, Minnesota, 1950, 8 p. 
(Transcript of a paper delivered at the Association for Computing 
Machinery Conference, March 28-29, 1950, Rutgers University), figs. 


An expository article (for a supplementary description of this information 
storage system see MTAC, v. 4, 1950, p. 31-39). 


5. CutHBert C. Hurp, The Selective Sequence Electronic Calculator of the 
International Business Machines Corporation, Carbide and Carbon Chem- 
icals Corporation, Engineering Development Division, Theoretical 
Analysis Department, Report No. K-134, Feb. 4, 1948, 6 pages, diag. 


An expository paper. 
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6. SHERMAN C. LOWELL, Digital Computer Research at the T. R. E., Malvern, 
Office of Naval Research, London, Technical Report ONRL—35-S50, 
Mar. 29, 1950, 3 p. 


The main design features of the Telecommunications Research Establish- 
ment computer are described. It is a general-purpose, high-speed, parallel, 
and binary machine. 


7. SHERMAN C. LOWELL, The Manchester University Digital Computer, 
Office of Naval Research, London, Technical Report ONRL—34-50, 
Mar. 29, 1950, 4 p. 


This is an expository article which describes a high-speed, serial, and 
binary machine and its applications and also a larger machine of similar 
design being built by Ferranti Brothers. The Ferranti computer will be 
superior from the engineering and maintenance standpoint. 


8. OFFICE OF NAVAL RESEARCH, Digital Computer Newsletter, v. 2, no. 4, 
December 1950, 5 p. 


The present status of the following digital computer projects is treated 
briefly in this number. 


1. NRL Computer (NAREC) 
2. Moore School Automatic Computer (MSAC) 
3. Naval Proving Ground Calculators 
4. SEAC (National Bureau of Standards Eastern Automatic Computer) 
5. Aberdeen Proving Ground Computers: 
Coupling of IBM Relay Calculators 


The ENIAC 
6. SWAC (National Bureau of Standards Western Automatic Com- 
puter) 
7. The Institute for Advanced Study-Computer 
8. Project Whirlwind 
9. The ACE Pilot Model, Teddington, England 
10. Data Handling and Conversion Equipment: 


Digital Reader 
Digital Data Recording System 
Analog Digital Converter: 


NEws 


American Statistical Association.—On Thursday December 28, 1950, at a meeting of 
the Association in Chicago, a round table discussion was held which had as its topic, ‘What 
can high-speed electronic computing equipment do for and to statistics?’’ The moderator 
was WILLIAM G. Mapow, University of Illinois. Mr. S. N. ALEXANDER, NBS, discussed the 
subject from the viewpoint of an electronic engineer while Mr. B. ScHREINER, A. C. Nielsen 
Company, talked about the problem as seen by an expert user. A discussion followed by H. 
C. Grieves, Bureau of the Census, and JoHN J. FINELLI, Metropolitan Life Insurance Co. 
Organizations joining in the meeting were the Psychometric Society, the Institute of Mathe- 
matical Statistics, and the Association for Computing Machinery. 
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Centre National de la Recherche Scientifique.—An international colloquium on “‘Mod- 
ern calculating machines and the human mind” was held at the CNRS on January 8 through 
13th at Paris. The program for the meeting was as follows: 


First Session, Recent Technical Progress in 
Large-Scale Calculating Machines 


“The Mark II, III, IV Machines” 


“A magnetic automatic calculating ma- 
chine” 

“Development of electronic digital com- 
puters at the National Bureau of Stand- 
ards” 

“Research activities of the National Phys- 
ical Laboratory” 

“The Institut Blaise Pascal Machine” 


“Mathematical machines in Sweden” 


“Research in progress at the University of 
Brussels” 

“Short-cut multiplication in a parallel-deci- 
mal automatic calculating machine” 


Second Session, A pplications of Mathematical 
and Scientific Problems to Large-Scale Cal- 
culating Machines 

“Numerical integration of the wave equa- 
tion” 

“Explanation of numerical methods of inte- 
gration of systems of linear partial differ- 
ential equations and the results obtained” 

“Different problems a machine can solve” 


“Les erreurs de chute dans les calculs sys- 
tematiques”’ 
“Operational experience with the EDSAC” 


Third Session, Large-Scale Calculating Ma- 
chinery and the Logic and Physiology of the 
Nervous System 

“A presentation of the homeostat”’ 


“Some new analogies between the structure 
of computers and the structure of the 
brain” 

“Large-scale calculating machines and the 
physiology of the nervous system”’ 

“Mechanical realization of models of the 
brain structure; presentation of artificial 
animals” 

“Work in Sweden” 


Louis DE BroG.tie, Institut Blaise Pascal, 
Secretary of Academy of Sciences, Chair- 
man 

H. H. Arken, Director of the Computation 
Laboratory, Harvard University 

A. D. Boots, Birkbeck College, University 
of London 

E. W. Cannon, National Bureau of Stand- 
ards, Washington, D. C. 


E. M. CoLesrook, National Physical Labo- 
ratory, Teddington, England 

L. CouFFIGNAL, Director of the Computer 
Laboratory, Institut Blaise Pascal 

Stic Exetorr, Chalmers Tekniska Hég- 
skola, Goteborg 

P. GerMAIN, Université Libre, Brussels 


E. J. Peruericx, S. H. HOoL.incpALce, 
Royal Aircraft Establishment, Farn- 
borough, Hants, England 

M. Caguot, Institut Blaise Pascal, Chair- 
man 


F. H. vaAN DEN DuNGEN, University of 
Brussels 
M. Piconr, University of Rome 


A. M. Uttiey, Telecommunications Re- 
search Establishment, Worcester, Eng- 
land 

A. VAN WIJNGAARDEN, Mathematische Cen- 
trum, Amsterdam 

M. V. WiikeEs, Director, University Mathe- 
matical Laboratory, Cambridge 

L. LapicguE, Institut Blaise Pascal, Chair- 
man 


W. R. Asusy, Department of Research, 
Barnhouse, England 

L. CouFFIGNAL, Institut Blaise Pascal, Di- 
rector, Laboratoire de Calcul Mécanique 


H. Gastaut, Faculté de Médecine de 
Marseille 

W. Grey WaLTER, Burden Neurological 
Institute 


G. KJELLBORG, Tekniska Hégskolan, Stock- 
holm 
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Discussion 


Discussion 
“Developments in automatic computing 
machinery in the Spanish school” 
“Présentation des appareils de Leonardo 
Torrés-Quevedo: 
1. le joueur d’échecs automatique 
2. le télékine, premier appareil con- 
struit pour le radio-guidage des 
bateaux 
3. les fusées logarithmiques de la ma- 
chine a résoudre les équations log- 
arithmiques” 


AUTOMATIC COMPUTING MACHINERY 


R. LorENTE DE No, Laboratories of the 
Rockefeller Institute 

W. McCuLtLocse 

P. Puic-Apam, Academy of Sciences, 
Madrid 

G. Torrés-Quevepo, Ingénieur des Ponts 
et Chaussées de Madrid 


Panel on Electron Tubes.—The Conference on Electron Tubes for Computers was 
held in Atlantic City on December 11th and 12th, 1950, under the joint sponsorship of the 
AIEE and the IRE in collaboration with the Panel. The program was as follows: 


Computer Experience with Electron Tubes 

“Review of AIEE conference on electron 
tubes for instrumentation and indus- 
trial use”’ 

“Experience with receiving-type vacuum 
tubes on the Whirlwind Computer 
Project” 

“Electron tube experience at IBM” 

“Performance of electron tubes in the 
ENIAC” 

Electron Tube Problems 

“Design and operation of tubes for long 
life” 

“The JETEC approach to the tube relia- 
bility problem” 

“Cathode inter-face impedance and its 
effects in aged vacuum tubes” 

“Cathode impedance and tube failure’”’ 

“6SN7WGE and 6AN5; mutual conduct- 
ance dispersion and the effect of low 
duty cycle operation on long life per- 
formance” 

“Open discussion in improved electron 
tubes for computers” 

Special Purpose Computer Tubes 

“A stable binary electrostatic storage 
system”’ 

“Recent experiences with the selective 
electrostatic memory tube” 

“The MIT storage tube” 

“The development of the Rogers Addi- 
tron” 

“A proposal for a binary adder tube util- 
izing beam deflection principles” 

“Special cold cathode discharge tube for 
counting and switching applications”’ 


Mrna REEs, ONR, Chairman 
W. R. Crark, Leeds and Northrup Co. 


E. S. Rica, MIT 


J. A. Goetz and A. W. Brooke, IBM 

Wricut E. Erton and Homer W. SPENCE, 
Aberdeen Proving Ground, Maryland 

S. N. ALEXANDER, NBS, Chairman 

J. O. McNatty, BTL 


J. R. Steen, Sylvania Electric Products, 
Inc. 


H. B. Frost, MIT 


L. S. NERGAARD, RCA 
I. Levy, Raytheon Manufacturing Co. 


A. L. SamuE., IBM, Chairman 
A. M. CLocston, BTL 


J. Raycuman, RCA 


P. Youtz, MIT 
T. VAN Dyk, Rogers Majestic, LTD. 


D. H. Griptey, NRL 


M. W. WaALLace and J. HEney, Federal 
Telecommunication Laboratories 





was 
f the 


NCE, 


ucts, 


sderal 


OTHER AIDS TO COMPUTATION 99 


“A decimal counting tube”’ 

“New improvements and applications in 
Remtron counter tube”’ 

Tube Manufacture and Crystal Diode Experi- 

ence 

“Design and manufacture of electron 
tubes for electronic computer service” 

“Problems in the manufacture of special 
tubes for computer usage”’ 

“Development of the 7AK7”’ 


“Some problems involved in the manu- 
facture of germanium diodes” 
“Experience with germanium diodes in the 
SEAC program” 
“Crystal diode life experience in the 
Whirlwind Computer circuits” 
Williams Type Storage 


“Tube experience in the SWAC” 

“The SEAC memory using Williams 
storage” 

“The selection of cathode-ray tubes for 
Williams storage”’ 

“Methods of testing cathode-ray tubes 
for service in Williams storage systems” 

“Theory of storage in cathode-ray tubes” 

“Progress report on the electron mecha- 
nism on the Williams storage tech- 
nique” 

“Space charge effects in Williams storage 
tube” 


T. R. Kouter, Philips Laboratories 
FRANK J. Cooke, Remington Rand Corp. 


J. G. Bratnerp, University of Pennsyl- 
vania, Chairman 


R. E. HicGs and H. E. Stumman, RCA 


R. L. McCormack, Raytheon Manufactur- 
ing Co. 

R. W. SLInKMAN, Sylvania Electric Prod- 
ucts, Inc. 

N. DeWotrr, General Electric Co. 


H. Wricut, NBS 

H. B. Frost, MIT 

J. H. BiGetow, Institute for Advanced 
Study, Chairman 

H. D. Huskey, NBS 

W. W. Davis, NBS 

J. H. PomMerene, Institute for Advanced 
Study 

D. FrreEpMAN, NBS 

J. Kates, University of Toronto 


A. W. Hott, NBS 


L. Bri_Ltourn, IBM 
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A New Differentiating Machine 


The first step in the design of a differentiating machine is the selection 
of a mechanical analogue for the derivative of a plotted function. The in- 
corporation of this mechanical differentiai ratio into a machine that will 
draw the derived curve of a plotted function will follow with the use of the 
proper linkages, gears, etc., that will transfer the relative value of this 
d 


. y — . . 
mechanical q to a writing pen for tracing the derived curve. 
x 


The most common type of differentiating machine that has been built 
uses a tangent line analogue based on the geometric concept of the deriva- 
tive. The differentiating machine built in 1904 by J. E. Murray! employs 
this tangent line analogue in the form of two dots on a celluloid plate. So 
long as the dots remain on the curve to be differentiated, the chord connect- 
ing them is approximately parallel to the tangent to the curve at the mid 
point between the dots, then by a system of connecting linkages, a writing 
pen is caused to draw the derived curve. 
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ARMIN ELMENDORF’s differentiating machine? employs a tangent line 
analogue in the form of a sharp edged wheel which is rolled along the curve 
to be differentiated. A pulley system of parallel linkages connects the wheel 
to a writing pen. Both Murray’s and Elmendorf’s machines can be adapted 
to use a ‘normal mirror’’ for finding the slope of a curve. A line perpendicular 
to a mirror held normal to a curve will be a tangent line. 

F. E. Myarp* based his machine on a different analogue from those 
already mentioned. His machine develops a differential ratio when a single 
point is made to trace a curve. This is done by an adding unit (a set of 
differential gears) and a crown wheel (integrating wheel) which together 

+n , . dy : a 
form the differential ratio .< . This ratio is plotted as a curve. 

No attempt is made in this article to evaluate the above mentioned 
machines as differentiating devices. 

A new differentiator has been designed and built by Cyrit P. ATKINSON 
(Fig. 1). It is a “tangent line analogue” type of differentiating machine 





Fic. 1. Atkinson’s differentiator. 


using a “‘tracer bar’’ to measure the slope of the tangent line. The relative 
value of the tangent of the angle of slope of the tracer bar is transferred to 
the writing pen by a pair of synchro-transmitters (selsyn motors) in such 
a way that the motor attached to the writing pen will cause it to describe 
the derived curve to a scale determined by the constants of the system. 

It had been hoped that the differentiator could be operated continuously 
in the differentiation of a curve and thus produce a continuous derivative. 
This plan had to be abandoned since an operator is unable to move the tracer 
bar with such finesse to cause the writing pen to draw a smooth curve. In 
the process of differentiation the machine is moved from interval to interval 
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across the original curve by means of the motor and drive shaft. At each 
interval a slope is determined by the tracer bar and a point is tapped by the 
writing pen, so that a series of points describes the derivative. A smooth 
curve drawn through these points is the derived curve. The units of the de- 
rived curve can be obtained from the scale factor of the machine. If the 
units of the abscissa of the original curve are equal to the distance between 
the centerline of the writing pen M of Fig. 1 and the center of rotation of the 
shaft of selsyn motor N, the units of the ordinates of the derived curve 
are equal to the ordinate units of the original curve divided by the units 
of the abscissa of the original curve. The abscissa units are the same for 
both curves. 

Various curves have been differentiated by the machine and the accuracy 
was determined to be from 2.5 to 5%. Fig. 1 shows the cosine curve under 
the writing pen of the differentiator. This cosine curve was produced by 
differentiating the sine curve shown under the tracer bar. Harmonic curves 
can be differentiated most easily when they are plotted with the units of 
the abscissa equal to distance between the centerline of the writing pen and 
the selsyn motor. When the slope of the tracer bar must exceed 45 degrees 
to be tangent to a curve as in the differentiation of a circular arc of 90 de- 
grees, a technique that is employed is to determine the normal to the curve 
by holding the tracer bar across the curve. By observing the reflection of the 
curve in the polished surface of the tracer bar and causing the reflection to 
be continuous with the original curve a very accurate normal can be ob- 
tained. Since a plot of the slope of the normal would be a curve of — = it 
is only necessary to take the negative reciprocal of this curve to obtain the 
derivative. 

Some of the defects of the earlier differentiators were kept in mind in the 
development of the present differentiator. For example, the tracer bar, the 
counterpart of which was obscured in both Murray’s and Elmendorf’s 
machines, is placed immediately before the operator and can be viewed from 
above and held tangent or normal to the curve without interference from 
other parts of the mechanism. Also, the writing pen, which is driven by a 
synchro-transmitter, can be placed in any position the designer desires, even 
in an adjoining room if need be. It was placed so the operator would have 
an unimpeded view of his work, both in the differentiation of the original 
curve and in the plotting of the derivative. 

The defects of the differentiator are (1) the point by point method of 
producing the derived curve; (2) the inaccuracies due to (a) the approxima- 
mation of the derivative by a chord, (b) selsyn motor inaccuracies, (c) 
mechanical inaccuracies. 

The machine could be improved by increased precision in the dimensions 
of the moving parts. The use of larger and more accurate selsyn motors 
would allow a derived curve to be drawn to a larger scale. The incorporation 
of photo tubes into the “tracer bar’ to replace the points which are used to 
trace the original curve would make the machine automatic and cause it 
to draw a continuous derived curve. 
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The comparison of the present machine with the machines mentioned 
earlier in this paper as to cost, performance, etc., and a more detailed account 
of the theory of differentiating machines awaits a more general article. 

C. P. ATKINSON 
hs von A. S. LEVENS 
Engineering Division 

University of California 
Berkeley 

1J. ErsKINE Murray, “A differentiating machine,” R. Soc. Edinburgh, Proc., v. 25, 
pt. 1, 1904, p. 277-280. (Reprinted in E. M. Horspurcu, Modern Methods and Instruments 
of Calculation. London and Edinburgh, 1914, p. 217-219.) 

2? ARMIN ELMENDORF, ‘Mechanical differentiation,” Franklin Institute, Jn., v. 185, 
am > 119-130. 

3. E. Myarp, ‘Nouvelles solutions de calcul grapho-mechanique—derivographes et 
planimeters,” Le Genie Civil, v. 104, 1934, p. 103-106. 

4MEYER zUR CAPELLEN, Mathematische Instrumente. Ann Arbor, 1947 (originally 
Leipzig, 1944). 

5 J. Lrpxa, Graphical and Mechanical Computations. New York, 1918. 

6 F. J. Murray, The Theory of Mathematical Machines. 2nd rev. ed., New York, 1948. 
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9. H. L. ANDREws, “‘Nomogram for G-M counter resolving time correc- 
tions,” Rev. Sci. Instruments, v. 21, 1950, p. 191. 


10. J. G. Bay ty, ‘‘An analog computer,” Rev. Sci. Instruments, v. 21, 1950, 
p. 228-231. 
d. 
A device is described for solving the equation = = — dx + f(t), where 
t is real time. The unknown function x is the rotation of a shaft driven by 


‘a 
a servomotor whose rate of rotation a regulated to equal a voltage 


d 
— rx + f(d). = is measured by a tachometer and the difference between it 


and — Ax + f(t) is used as an input for a servo-system of the type with a 
“memory circuit” described by W1LLIaMs & UTTLEY.! A helical potentiom- 
eter on the output shaft yields the voltage Ax and f(t) is introduced by a 
manual follower. Various applications to atomic pile calculations are 
described. 

F. J. M. 


1F, C. Wittiams & A. M. UtTTLey, “The velodyne,” Inst. Elec. Engrs., Jn., v. 93, 
part IIIA, 1946, p. 1256-1274. 


11. A. BEtsErR, “A slide rule for nuclear emulsion calculations,’’ Rev. Sci. 
Instruments, v. 21, 1950, p. 933-934. 


12. L. M. Haupt, “Solution of simultaneous equations through use of the 
A.C. network calculator,” Rev. Sci. Instruments, v. 21, 1950, p. 683-686. 


The system of equations to be solved are realized by means of Kirchhoff’s 
laws. The coefficients are represented by impedances either directly or by 
reflection through one-to-one transformers. This last device is used to remove 
the restrictions which normally apply to real transfer impedances. 

F. J. M. 
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13. W. A. McCoo1, ‘‘Frequency analysis by electronic analog methods,” 
NRL Report no. 3724, Aug. 21, 1950. 


The complex Fourier transform F(w) of the real function f(r) is ex- 
pressed as 


Flw) = {- flrye-*"dr. 


F(w) may be evaluated on analog equipment by a solution of the linear 
integrodifferential equation 


ytu fxt=s0; 0) = f yt =. 


From a solution of a similar system of equations one may obtain the real 
inverse Fourier transform of F(w). An illustrative example and an analog 
machine program are given. The author has reported several omissions and 
errata that do not affect the results. 

PauL Brock 


Reeves Instrument Corp. 
New York, N. Y. 


14. RicHarD MCFEE, “‘A trigonometric computer with electrocardiographic 
applications,” Rev. Sci. Instruments, v. 21, 1950, p. 420-426. 


The author points out the usefulness of a double electronic resolver which 
will give the three spherical coordinates of a vector from its three cartesian 
coordinates expressed as voltages. Such a spherical resolver is readily ob- 
tained by combining two ordinary resolvers of the two dimensional type. 
An electronic resolver which yields (e;? + e,”)! and the arctan é2/e; is de- 
scribed. The input voltages e, and e2 are used to modulate two sinusoidal 
voltages of the same frequency, 90° out of phase. These modulated voltages 
are added and the amplitude is detected to yield (e,* + e,*)!. The phase of 
this sum voltage is also detected by a pulse technique. A pulse is emitted 
when the reference voltage passes through zero and also one when the sum 
voltage passes through zero. The difference in time between these is the 
required angle arctan é2/e;. Special provision is made to yield a zero angle 
for zero amplitude and for obtaining the base line for amplitude measure- 
ments. A proposal is also made for displaying the results in the three dimen- 
sional case. 


F. J. M. 


15. F. G. FENDER, “‘Aerocom. . . . An analog computer,”’ Northwestern 
Engineer, v. 9, no. 1, Mar. 1950, p. 10-12, 32. 


This article describes a large analog computer located in the Aerial 
Measurements Laboratory at Northwestern Technological Institute. This 
device “‘now consists of the following types of equipment: inductors, resis- 
tors, and capacitors; special amplifiers; special synchronous switches and 
relays to provide repetition of the transients and to re-establish initial 
conditions between such repetitions; various function generators; single and 
dual gun cathode ray oscilloscopes to provide visual representation of out- 
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puts together with cameras for making permanent records; and various 
standard electrical devices for calibration and testing. Circuit connections 
are made by busses with selector switches and by plugs and jacks. . . .” 
The cost was about $250,000. 

F. J. M. 


16. R. L. Garwin, “A differential analyzer for the Schrédinger equation,” 
Rev. Sci. Instruments, v. 21, 1950, p. 411-416. 


The equation y’”’ — [E — V(t) ]¢ = 0 is solved by means of feed back 
amplifier integrators. The potential V(t) is obtained by a manual follower. 
The design considerations for a simple and relatively inexpensive instrument 
are given at some length and also its use for the solution of the appropriate 
characteristic value problems. 

F. J. M. 


17. B.O. MARSHALL, JR., ‘““The electronic isograph for roots of polynomials,” 
Jn. Appl. Phys., v. 21, 1950, p. 307-312. 


A polynomial f(z) = >> a;z‘ with real a; can be represented as a function 
j=l 
of R and @ when z = Rexp #, i.e., f(x) = >> a;R/ cos j@ + i ¥ a;R’ sin 76. 
7=1 j=1 
If R is entered manually as a shaft rotation, it is possible to obtain d.c. 
voltages corresponding to the quantities a;R’ by a suitable arrangement of 
potentiometers. These d.c. potentials are applied to commutators to give a 
square wave voltage of frequency 607 per second. All but the fundamental of 
this voltage is filtered out and added to yield the real and imaginary parts 
of f(z). These are used as the horizontal and vertical deflection voltages of 
an oscilloscope. R is varied until a zero is observed. The instrument handles 
polynomials of the tenth degree and the output voltages have three figure 
accuracy under suitable scaling. 
F. J. M. 


18. FRANCGOIS-HENRI RAYMOND, “Sur un type général de machines mathé- 
matiques algébriques,”’ Ann. Telecommun., v. 5, 1950, p. 2-19. 


The author considers three particular mathematical machines of the 
continuous type. One is a machine for solving systems of linear algebraic 
equations and associated problems (e.g., inversion of a matrix, characteristic 
values of a matrix). This will be called the ‘‘algebraic machine.’’ The other 
two machines are for solving systems of differential equations. One machine 
is capable of handling larger systems than the other. However, the theory 
of both is the same and either machine will be referred to as a ‘‘differential 
analyzer.’’ After some introductory remarks the author briefly describes the 
operation and applications of the machines. The next portion of the paper, 
which is of greatest interest to mathematicians, is a discussion of the sta- 
bility and precision of the devices. These phases will be elucidated below. 

Consider the system of linear equations (in matrix form) Ax = b where 
A = |la,;|| is a square matrix and x and 6} are column vectors. It is assumed 
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that by a preliminary manipulation all the coefficients of A have been re- 
duced (by a scale factor) such that |a,;| = 1. It is also assumed that to each 
coefficient a;; there is attached an error a;; whose upper bound is a certain 
percentage of a;;. Similarly for the b; (components of the vector 5). Hence 
the machine actually solves the equation (A + a)x’ = 6”. The problem is 
to estimate |x,’ — x;,|, that is the difference between the observed solutions 
x,’ and the actual solutions x;. By using the REDHEFFER formula! the author 
obtains the following result: 


As eD-“(L + 2'Xy)(T/(n — 1))*? 


where A = max|x,’ — x;|, € = percentage of the coefficient a;; which forms 


the error aj;; i.e., aj; = cai; and |b,’ — b;| = €|b;|, D = the modulus of the 
determinant of A, T = the trace of AA, (A; = transpose of A), L = the 


length of the vector 6 = [> b?7]!, Xu = max|x,’|, » = number of equa- 
i=1 é 


tions in the system. 

The quantity A measures the percentage error. However it appears to 
the reviewer that the scaling operations are not considered in sufficient 
detail; for while the percentage error may be small, the actual error may be 
large if a large scale factor has been used. Also the justification for the as- 
sumption that a;; = ¢a;; is not obvious. Noise is not considered. 

The algebraic machine is essentially of the GOLDBERG-BROWN type’; and 
the stability of the system depends on the characteristic roots of the matrix 
A. In writing the equilibrium equations for the machine, the following 
matrix form is obtained: 


(1) [A — ((n + 1)/G)I x —b =0 
where G(p) is the gain characteristic (in operational form) of the amplifiers 
used, (p = d/dt is the Heaviside operator). If u,, » = 1, 2, ---, m are the 


characteristic roots of A, then uw, = (m + 1)/G(p). If the amplifiers are so 
constructed that RG(p) < 0 when RP > 0 then Rp > O if Ryu, > 0 and the 
machine is stable. [R = “real part of.’’] 

In the differential analyzer the equations to be solved are (in matrix form) 


(2) g(t) + [Ap* + Bp? + Cp + D]Y = 0 


where A, B, C, and D are square matrices with constant coefficients, g(t) 
and Y are column vectors and » = d/dt. In a manner similar to that used in 
deriving (1), 


{[A — ((n + 1)/G)I] + Ba + Co? + Da®}V+2=0 
is obtained. [In the ideal case, G = © (infinite gain amplifier) and 


a(p) = —p-.] As in the algebraic case, everything depends on the charac- 
teristic roots of the matrix 


(3) A + Ba + Ca’? + Da’. 


In the ideal case, the characteristic frequencies are the roots of the determi- 
nantal equation 
det|A — Bp — Cp — Dp-| = 0. 
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Call them ~;. The characteristic roots of the matrix (3) will differ from the 
px by small amounts, 5p, which the author computes later (using the first 
few terms of a Taylor series expansion). However this matrix (3) may have 
extraneous roots, g,, vy = 1,2,-:-. 

Under the assumptions |G(g)| <|G(0)| and |a(p)| is small (~10-*) it 
can be shown that the g, depend only on the matrix A and not on B. The 
author states that if this is not the case his results are not valid and that a 
general discussion appears difficult and is outside the scope of the present 
study. 

Making certain (practical) assumptions regarding the physical system, 
for example, assuming G(p) = — Go/(1 + pT) where Gp = |G(0)| and also 
assuming (for the present) that A is symmetric, the author obtains the results 
that two sufficient conditions for the correct functioning of the differential 
analyzer are (a) |p| < 2rAfo, (b) |qm| > 2rAfo where |gn| = min|g,| and 


Af is the band-width of the amplifier (27Af,>7) = 1). Using the Redheffer 
formula again, it is also shown that 


Qu = — [1 + Go/(m + 1)||A||((m — 1)/T)"*]T 


where ||A|| = |det A| and T is the trace of A. 

The author remarks that (a) can always be satisfied by an appropriate 
change of scale of the independent variable in the original system of equa- 
tions. Since Go is always large (~10*) in a well designed amplifier, he con- 
cludes that the critical parameter, from the point of view of stability and 
precision, is ||A||/7*—. Errors in the matrices A, B, C, D are not considered 
(neither steady state nor random). 

Certain results are obtained if the condition that the matrix A be sym- 
metric is relaxed. 

K. S. MILLER 


New York University 
New York 


1R. REDHEFFER, “Errors in simultaneous linear equations,” Quart. Appl. Math., v. 6, 
1948, "4 342-343. 
2 


A. GotpBErRG & G. W. Brown, “An electronic simultaneous equation solver,” 
Jn. Appl. Physics, v. 19, 1948, p. 339-345 [MTAC, v. 3, p. 329-330]. 


19. A. H. Scott, “An instrument for mechanically differentiating curves,” 
Rev. Sci. Instruments, v. 21, 1950, p. 397-398. 


A Corap! intergraph is modified by replacing the pen by a stand hold- 
ing a slide rule glass slide. The method of operation is carefully described in 
this note. An accuracy of two percent in the permissible range of slopes 
was obtained. 


F. J. M. 


20. H. R. SEIWELL, “‘A new mechanical autocorrelator,”’ Rev. Sci. Instru- 
ments, v. 21, 1950, p. 481-484. 
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by means of two ball cage integrators. f,; and f, are entered by manual 
followers from graphs. f2 is utilized as the displacement of the cage of the 
first integrator whose disk is driven at a constant rate. The output of the 
first integrator drives the disk of the second and the cage displacement of 
the second integrator is proportional to f,(#). The output of the second 
integrator is registered on a counter and is the required answer when the 
interval of integration has been covered. The accuracy of the device depends 
on the slope of the functions f. 


F. J. M. 


21. H. Sarmizu, P. J. Etsey & D. McLacaian Jr., “A machine for syn- 
thesizing two-dimensional Fourier series in the determination of crystal 
structures,” Rev. Sci. Instruments, v. 21, 1950, p. 779-783. 


The stators of a pair of selsyns are connected so that a voltage across 
the rotor of the first induces a voltage in the rotor of the second proportional 
to the cosine of the sum of the angles of rotation of each rotor from a refer- 
ence position. Thus if the first is rotated an amount mx and the second an 
amount my a term Gam cos (nx + my) can be represented. These voltages can 
be immediately added to produce F(x, y) = >> dam cos (nx + my). The co- 
ordinates x and y are entered as the rotation of two shafts and the rotations 
nx and my are taken from these by a suitable arrangement of gears. In the 
present device m and m each have a range of 8 but it is planned to extend 
this range to 16. 


F. J. M. 


22. F.C. SNOwpEN & H. T. Pace, “‘An electronic circuit which extracts 
antilogarithms directly,” Rev. Sci. Instruments, v. 21, 1950, p. 179-181. 


This circuit is based on the use of a 6SK7-GT as an inverted triode. The 
grid current is then proportional to the antilogarithm of the applied plate 
voltage for a range of 15 volts and for a range of grid current from .33 to 
1 ma. 


F. J. M. 


23. Rospert R. Rem & Du Ray E. Strompack, “Mechanical computing 
mechanisms,” Product Engineering, v. 20, 1949, August no. 8, p. 131-135, 
Sept. no. 9, p. 119-123, Oct. no. 10, p. 126-130, Nov. no. 11, p. 121-124. 


In the first article of this series of four, specifications for the input and 
output of mechanical computers are described. Three types of errors are 
described: Class A errors due to inaccuracies in components. Class B errors 
due to mathematical incompleteness in the setup and Class C or personal 
errors. It seems to be desirable to accept certain Class B errors which are 
known in order to minimize Class A and Class C errors. Designs questions 
in relation to scaling are discussed. In the second article, cams, resolvers, 
trigometric function generators and differentials are discussed relative to 
available scales and error characteristics, in the above classification. The 
third article deals with multiplication, division, integration differentiation 
and link mechanisms. The last article describes a number of applications. 

F. J. M. 
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24. HENRY WALLMAN, “An electronic integral transform computer and the 
practical solution of integral equations,’ Franklin Inst., Jn., v. 250, 
1950, p. 45-61. 

A proposed device is described for presenting 


¢ K(x, t)f()dt 


in the form of a graph of a function of x on a cathode ray tube. K(x, t) is to 
be obtained by scanning a photographic plate whose opacity corresponds to 
the value of K at the point x, ¢. The multiplication of K and f and the value 
of f(t) itself are to be obtained by using the components due to MACNEE.' 
The author shows how the persistence of the image on the screen of the 
cathode ray tube can be utilized to construct an iterate of the transform and 
describes also how non-linear transforms in the form 


'b 
f Kee, ore, soya 
can be obtained. 

The paper also lists a variety of applications of such a device. The case 
in which K(x, t) = cos xt yields the impulse response of an electrical net- 
work. Such a device would readily yield the coefficients of the orthogonal 
expansion of a function, the Hilbert transform and the convolution integral. 

Another set of applications is concerned with the solution of the integral 
equations. Special cases treated include simultaneous linear algebraic equa- 
tions, the Volterra equation, the use of Liouville-Neumann series, Fred- 
holm’s integral equation of the second kind, the Dirichlet problem for a 
plane potential and a non-linear problem for a pendulum. 

F. J. M. 


1A. B. Macnee, “‘A high speed electronic differential analyzer,” I.R.E. Proc., v. 37, 
1948, p. 1315-1324 [MTAC, v. 4, p. 119-120]. 


NOTES 


124. LESLIE JOHN CoMRIE (1893-1950).—This great table maker and 
pioneer in the art of mechanical computation was born in New Zealand in 
1893. He received his early training and a M.A. degree at the University 
of New Zealand. He saw active service during the first world war with the 
New Zealand Expeditionary Forces and after the armistice went to Univer- 
sity College, London and Cambridge University, where he received his 
Ph.D. in Astronomy in 1923. After 3 years teaching in the United States at 
Swarthmore College and Northwestern University he returned to England and 
the Royal Greenwich Observatory as Deputy Superintendent of H.M. 
Nautical Almanac Office. He became Superintendent in 1930 and held that 
post until 1936. Here he introduced modern computing methods which did 
much to increase the efficiency and productivity of the office. It is this work 
which brought out his genius for the organization and keen analysis of com- 
puting and table preparation for which he later became so famous. He also 
served brilliantly as secretary of the BAASMTC during 1929-36 and was 
much concerned with the production of the committee’s first six volumes. 
In 1937 he left the Observatory to devote his entire energy to the develop- 
ment of the Scientific Computing Service, the first enterprise of its kind. 
The history of this organization is one of lasting achievement and pioneering 
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effort. Its contributions to the war effort in the 1940’s was the source of 
much justifiable pride to Comrie. He served as cooperative editor of MTAC 
from 1944-49, and in 1950 was elected Fellow of the Royal Society. His 
death, Dec. 11, 1950, came as a relief from a lingering affliction. 

The reader is referred to Mathematical Table Makers' for a complete 
account of Comrie’s writings and tables. He is perhaps known most widely 
as the editor of the third (1930) and later editions of Barlow’s Tables. His most 
recent work was the preparation of the monumental two volume edition of 
Chambers’s Tables. He delighted in the minute analysis of a computing prob- 
lem with respect to a given machine and took great pleasure in finding scien- 
tific uses for a particular feature of a machine which had been intended by 
the manufacturer for some trivial commercial application. This exploitation 
of the commercially available equipment with its built-in, mass-production 
precision and service was always uppermost in his mind. He had no great 
enthusiasm for the unreliable specially made computing device. | can recall 
many friendly arguments with him beginning in 1932 on these two opposing 
points of view, my side of the argument being, in those days, rather difficult 
to maintain. His knowledge of numerical analysis was profound and at the 
same time severely practical. He will be remembered for his ‘‘throw back”’ 
method of modified differences. On his last visit to America in 1946 he was 
impressed by what he saw of new computing techniques, but also somewhat 
dismayed to see to what uses they were put. At one point of the inspection 
he took me aside to remark: ‘‘These people have a terrific amount to learn 
about computing.”” This remark is as applicable today as it was when he 
must have made it (to himself) as a young man at Greenwich. It seems to 
have been the keynote of this crusading calculator. 

D. H. L. 


'R. C. ARCHIBALD, Mathematical Table Makers (The Scripta Mathematica Series no. 3), 
New York, 1948 [MTAC, v. 3, p. 143]. In this work will be found two portraits of Comrie. 


125. THE d*? Test oF RANpom Dicits.—The testing of digits for local 
randomness, when the digits are to be used in a Monte Carlo method, seems 
to require a different type of test from the four proposed by KENDALL and 
BABINGTON-SMITH.! Kendall’s tests (frequency, serial, gap, poker) apply to 
random digits as used normally in random sampling, etc. 

In a Monte Carlo method, random digits are used to select a random 
point in the unit square, the digits thus representing the coordinates of a 
point between (0, 0) and (1, 1). 

For two such points, the probability that the square of the distance 
between them be less than a? is given by? 


4 
P = re? —- o. 
for a® = 0.0, 0.1, 0.2, ---, 0.9 


8 4 
P = 3+ (x — 2)a® + 4(a? — 1)! + 3 (a? — 1)! - 7 — 4a? sec“ a 
for a* = 1.0, 1.1, 1.2, ---, 2.0. 


P is, of course, a continuous function. Discrete values of a? are selected 
for convenience in computation. 
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These probabilities are, in order: 


oe } a P 
0.1 .234832 1.1 .985703 
0.2 .409805 1.2 .992048 
0.3 .549300 4.3 .995788 
0.4 .662018 1.4 .997926 
0.5 .752987 1.5 .999080 
0.6 .825601 1.6 999652 
0.7 .882349 a7 .999898 
0.8 .925163 1.8 .999982 
0.9 -955593 1.9 .999999 
1.0 -.974926 2.0 1.000000 


In applying the test, sets of 3 digits are sufficient to represent one co- 
ordinate; i.e., 12 digits are selected to give both coordinates of two points. 

In a limited test, no significant difference was found when using two-digit 
coordinates; in other words, the third digit has only a minor effect on the 
square of the distance. The square of the distance is computed and punched 
to one decimal place. The distribution of the results is compared to the 
theoretical by the Chi-squared test. The calculation can be readily carried 
out on any of the IBM calculators (602, 602A, 604). The distribution of the 
results can be made in one run on a tabulator equipped with two digit 
selectors. 

This test was carried out on the random digit table of the University of 
Wisconsin Computing Service (10,000 cards, each bearing 40 digits), making 
six calculations of d? on each card, selecting the columns to be read at 
random. Results were tabulated for each 1000 cards (6000 d?’s); the x? 
analysis giving: 


Thousand x’ p 
1 17.085 .26 
2 11.847 62 
3 18.585 19 
4 7.529 91 
5 7.762 91 
6 11.340 .66 
7 14.162 44 
8 19.088 17 
9 8.370 .86 

10 9.747 81 


using 14 degrees of freedom, since the last 6 classes were lumped together. 

Wiring diagrams for the d? calculation on the 602A calculating punch 
and the distribution of results on the 416 or 405 tabulator are available 
on request. 


. ; “RED GRUENBERGE 
Computing Service FRED GRUENBERGER 
University of Wisconsin 
Madison, Wisconsin 


A. M. MARK 


Southern Illinois University 
Carbondale, Illinois 
1M. G. KENDALL & B. BaBincton-SmitTH, R. Stat. Soc., Jn., v. 101, 1938, p. 157, and 


Supplement, v. 6, 1939, p. 51. 
2? BENJAMIN WILLIAMSON, Integral Calculus. 6th ed., London, 1891, p. 390. 
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126. ON THE CALCULATION OF THE SQUARE Root By AUTOMATIC Com- 
PUTING MACHINES.—The calculation of the square root of a number on an 
automatic computing machine not equipped with an internally controlled 
square root order requires the application of some standard sub-program. 
It is generally convenient to employ a well-known iterative procedure based 
upon the formula 


(1) Ro) = LR + N/R], 


where R™ is the nth approximation to the square root R = N!. The efficient 
application of this procedure requires a systematic method of selecting a 
zeroth approximation R® to the square root R, and the number of iterations 
required to calculate R to a specified number of significant figures depends 
largely upon the accuracy of the zeroth approximation. 

Let N be a decimal number in the range 1 < N < 100. Define — by 


R=é£+N/10, 
and let the zeroth approximation R® be given by 
R® = + N/10, 
where £ is an average value of £. If the number JN is distributed uniformly on 
a logarithmic scale,' we may take 
of 100 
— = (log 100) f (N! — N/10)d log N = 1.759 ~ 2, 
1 
the one-digit approximation — = 2 being sufficiently accurate. Therefore, 
we let 
(2) R® = 2 + N/10. 


With (1), the first approximation R® may be written 


1fN 10N 
RY =5/ 5 +2+ 


10 N + 20 
This approximation differs from the root R by an amount A = R — R, 
where 
et 
7 20 © R? + 20 


In the range 1 < R <¢ 10, A is positive. It is equal to 0.25 at R = 1, de- 
creases to a minimum of 0.038 at R = 1.957 --- and is equal to 0.166 --- at 
R = 10. Therefore, the choice, (2), of a zeroth approximation leads to a 
first approximation with at least one correct significant figure. The calcula- 
tion of the square root by the application of (1) correct to eight significant 
figures requires a maximum of five iterations and to sixteen significant figures 
a maximum of six iterations. 

If N* is any positive number such that N* = 10*%N,1 < N < 100, the 
zeroth approximation R®* to the root R* = (N*)! is 


R* = 10*(2 + N/10). 
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This method has been applied in this laboratory to the formulation of a 
square-root order on control panels for 10-digit arithmetic utilizing the IBM 
Card-Controlled Electronic Calculator. 

ROBERT W. SMITH, JR. 
Stuart R. BRINKLEY, Jr. 
Explosives and Physical 
Sciences Division 
U. S. Bureau of Mines 
Pittsburgh, Pa. 

1 It is reasonable to assume that dimensional numbers will be uniformly distributed on 

a logarithmic scale if the choice of dimensional units is random. However, this point is 


unimportant, since this result is comparatively insensitive to the nature of the averaging 
process. 


QUERIES 


37. THE SQUARE Root METHOD FOR LINEAR EQuaTIONs.—In a letter 
dated 7 Feb. 1950 Mr. H. F. Ratnsrorp of Colonial Surveys, Bushy Park, 
Teddington, England, commenting on the article entitled ‘‘The square root 
method for solving simultaneous linear equations’ by J. LADERMAN in 
MTAC, v. 3, p. 13-16, points out that this method was not “‘probably first 
discovered by Banachiewicz in 1938’ but goes back at least to CHOLESKY 
whose treatment of the problem was described by BENoit' in 1924. Can any 
reader supply an earlier reference to this method? 

1 Benoit, “Note sur une méthode de résolution des équations normales provenant de 
l’application de la méthode des moindres carrés 4 un systéme d’équations lineaires en nombre 
inférieur 4 celui des inconnues. Application de la méthode a la résolution d’un systéme define 
d’equations linéaires,” International Geodetic and Geophysical Union, Association of 
Geodesy, Bullétin Géodésique, no. 2, 1924, p. 67-77. An English translation of this article has 
been kindly supplied by Mr. Rainsford and is available in the UMT FILE. 


QUERIES—REPLIES 


47. RusSIAN BESSEL FUNCTION TABLEs (Q 25, v. 3, p. 66).—The volume 
referred to in this Query, namely: Tablifsy Znachenit Funkisit Besselia ot 
Mnimogo Argumenta, was not published until 1950. [It will be reviewed in 
the next issue of MTAC.] 

R. C. ARCHIBALD 
Brown University 
Providence, R. I. 
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CORRIGENDA 113 
; “a CORRIGENDA 
V. 4, p. 130 equation (2) for superscripts p + 1 and p + 2 read p + 2 and p + 4. 
V. 4, p. 151, 1. —4, for 23 read 33. 
R. V. 4, p. 205 1. 2, for by. read by:. 
V.4, p. 244, 1. 1 for on read an. 
V. 4, p. 255-259. The following omissions and corrections in the Name Index were kindly 
supplied by S. A. Jorre: 
Insert 
_—. Aiken H. H., 228 Lipkis R., 108 
wr fon Archibald R. C., 124 McPherson J. L., 117 
Bernoulli J., 189, 191, 208-212 Mobius A. F., 84 
Bragg W. L., 177 NBSCL, 161, 163 
Buchner J. P., 191 Nye J. F., 177 
Clausberg C., 198 Pell J., 189 
Cossar J., 179 Rao C. R., 209 
letter en 4 <4 a esc pt 7S 
Park, erdélyi A., Baie Scherberg M. G., 173 
3 t Heisenberg W., 151 Sharp A., 198 
c ~~ Herget P., 132 Sherwin H., 198 
AN in Hinrich J. C., 192 Smith C. S., 16 
y first Hobson E. W., 185 Stibitz G. R., 114 
JESKY Huskey H. D., 108 Stone A. H., 98 
n any Jahnke E., 152 Thomas T. Y., 99 
Kolmogoroff A., 118, 127 Tricomi F. G., 217 
wit Kiissner H. G., 153 Vega G., 199 
ant de ‘ aha , 
alien Lévy P., 204 ’ Whitfield J. W., 210 
. define Liebmann H., 75 Wolfram I., 191 
ra of Liénard A., 94 
cle has 
Delete 
Késsner, Sherberg, Smith C. B. 
For Read For Read 
Erdélyi Erdélyi Stieljes Stieltjes 
euIS Posch Pésch Vantil Van Til 
lia ot Prevost Prévost Weiner J. P. Weiner J. R. 
red in Spagenberg Spangenberg Zagor H. E. Zagor H. I. 
Shapley H., 26 Shapley H., 25 
LD 
Interchange 
Dennis and Dempsey, Harrison and Harris, Peterson and Peters, Zuckerman and Zucker 














